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Abstract. An abstract framework for the theory of statistical solutions is devel¬ 
oped for general evolution equations, extending the theory initially developed for 
the three-dimensional incompressible Navier-Stokes equations. The motivation for 
this concept is to model the evolution of uncertainties on the initial conditions for 
systems which have global solutions that are not known to be unique. Both con¬ 
cepts of statistical solution in trajectory space and in phase space are given, and 
the corresponding results of existence of statistical solution for the associated initial 
value problems are proved. The wide applicability of the theory is illustrated with 
the very incompressible Navier-Stokes equations, a reaction-diffusion equation, and 
a nonlinear wave equation, all displaying the property of global existence of weak 
solutions without a known result of global uniqueness. 


1. Introduction 

The concept of statistical solution was introduced for the study of turbulence in 
incompressible Newtonian fluid flows. In a turbulent flow, most relevant physical 
quantities (e.g. velocity, kinetic energy, energy dissipation) display a wild variation 
in space and time, while displaying a more orderly behavior when averaged in space 
or time (see e.g. [68, 5, 42, 54, 38, 51]). This behavior appears, in fact, for different 
realizations of the flow, with somehow “universal” properties, so that one is led to 
consider averages with respect to an ensemble of flows, in an attempt to capture 
common properties of turbulent flows. 

One can work with ensemble average in a formal sense, without worrying about the 
regularity of the solutions of the system, or one can be more strict and work with some 
notion of weak solution of the system for which existence results for the corresponding 
initial value problem are available. The statistical solution were introduced exactly 
with this later purpose in mind: they have been defined to model, in a rigorous 
way, the evolution of ensembles of weak solutions of the incompressible Navier-Stokes 
equations, as a foundation for a rigorous treatment of turbulent flows. 

Since its inception, the theory of statistical solutions for the incompressible Navier- 
Stokes equations has been the basis for a growing number of rigorous results for 
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turbulent flows (e.g. [29, 22, 8, 30, 39, 32, 31, 60, 63]). The concept of statistical 
solution has also been successfully adapted to a number of other models, particularly 
fluid flow models, but also other types of nonlinear partial differential equations (e.g. 
[14, 4, 77, 44, 67, 19, 74, 58, 17, 16, 23, 47, 56, 24, 46, 61, 13]). 

In fact, the notion of ensemble average is relevant for any evolution system dis¬ 
playing a complicated dynamics, in which uncertainties in the initial condition are of 
crucial concern. The concept of ensemble averages is directly related to the evolu¬ 
tion of a probability distribution of initial conditions. In a well-posed system, with 
a well-defined semigroup {S(t)} t > 0 , the evolution {/r t } t >o of the probability distribu¬ 
tion of the state of the system at each time t is just the transport, or push-forward, 
Ht = of the initial probability measure /x 0 , by the semigroup (more precisely, 

Ht(E) = (t) -1 E), for any Borel subset E of the phase space). The difficulty is 

to extend this definition to obtain the distributions /q for systems in which {S(t)}t >o 
might not be defined. This was the aim of the concept of statistical solution in the par¬ 
ticular and fundamental case of the three-dimensional incompressible Navier-Stokes 
equations. 

Two main definitions of statistical solutions have been introduced in the 1970’s. 
First, Foias ([27, 28]), in works stemmed from discussions with Prodi (see e.g. [34]), 
introduced the concept of statistical solution in phase space, consisting of a family 
of measures on the phase space of the Navier-Stokes system, parametrized by the 
time variable, and representing the evolution of the probability distribution of the 
state of the system. Then, Vishik and Fursikov ([75, 76]) introduced the notion of a 
space-time statistical solution, which is that of a single measure defined on the space 
of trajectories of the system, hence encompassing both space and time variables at 
the same time. 

Still in the mid to late 1970’s, it is worth mentioning the work by Ladyzhenskaya 
and Vershik ([49]), presenting a different proof of existence of statistical solution (in 
phase space) using a representation theorem by Castaing ([21]), for the measurability 
of multivalued solution maps, to overcome the fact that the solution operator is not 
a continuous, single-valued map. There is also the work by Arsen’ev ([2]) using a 
measurable selection argument to construct a statistical solution and, in fact, already 
introducing a notion of space-time statistical solution, which at first had some restric¬ 
tions on the initial measure, but that were eventually relaxed ([3]). Another proof 
of existence of space-time statistical solutions, based on nonstandard analysis, was 
given in the early 1990’s, by Capinski and Cutland (see e.g. [15]). 

More recently, inspired by the definition given in [75, 76] and using a different 
approach for the construction of statistical solutions, by approximating the initial 
measure by convex combinations of Dirac delta functions, given in [33], Foias, Rosa 
and Temam [35] considered the idea of space-time statistical solution of Vishik and 
Fursikov but with slightly different hypotheses that make it more amenable to anal¬ 
ysis. Projecting that modified space-time statistical solution to the phase space at 
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each time yields a family of measures which is a particular type of statistical solution 
in phase space, hence bridging the two notions given earlier. 

This work [35] inspired us to look for a more general abstract framework for the 
definition of statistical solutions for various types of systems and with reasonable 
conditions guaranteeing the existence of statistical solutions for the associated initial- 
value problem. The constructions developed in the previous works were very much 
based on the structure of the equations, but nevertheless we hoped to find a framework 
general enough to be applied to different types of systems. 

The first result was published in [12] and was sufficiently general to apply naturally 
to an ample class of problems. However, there is one condition in [12] which has a 
more involved formulation and made us feel that there should be a simpler way to 
attack the problem. Moreover, [12] only addressed the notion of space-time statistical 
solution, which we call here a trajectory statistical solution or a statistical solution 
in trajectory space. 

In our current work, we keep the same framework of [12] for the trajectory statistical 
solutions but with a simpler set of hypotheses that avoid the more involved condition 
in our previous work. The proof is also simpler with this new set of hypotheses. 
Furthermore, we address not only statistical solutions in trajectory space but also 
in phase space, presenting a general setting for the definition of those two types 
of solutions and the corresponding theorems of existence for the associated initial 
value problems. The proof also significantly simplifies the known proofs for concrete 
examples of equations, particularly the Navier-Stokes equations. 

We should mention that the concept of statistical solution had been previously 
extended to a few frameworks encompassing some particular classes of differential 
equations (e.g. [43, 40, 7, 41, 66]), but none of them nearly as general as done here. 

The hypotheses needed for our existence results are natural and not difficult to 
verify, relying essentially on properties of the set of individual solutions of the system. 
In fact, the key ingredient in our approach is to look for topological and measure- 
theoretic properties of the set of solutions instead of looking at the structure of the 
equation. 

The abstract framework that we construct starts with a Hausdorff space X, an 
interval / C R, the space of continuous paths X = Ci oc (/, X) endowed with the 
compact-open topology, and a subset U of X. There is no system of equations or so¬ 
lution operator at this abstract level, allowing for a wide range of applications, even 
for evolution problems not arising directly from differential equations. Only in the 
applications is that those objects will be realized, with the space X being interpreted 
as the phase space of the system; the interval I, as a time interval for the evolution 
system; the space X, as a space-time function space in which the solutions, or tra¬ 
jectories, of the system are included; and U , as a subset of solutions, or trajectories, 
of the system. In the case of partial differential equations, the set U is usually the 
whole set of weak solutions of the problem or some particular subset of solutions 
(e.g. Leray-Hopf weak solutions or suitable weak solutions for the three-dimensional 
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Navier-Stokes equations, viscosity solutions for evolutionary Hamilton-Jacobi equa¬ 
tions). 

In this abstract setting, a ^/-trajectory statistical solution or simply a trajectory 
statistical solution is a tight Borel probability measure p carried by a Borel subset of 
X which is included in the set U (see Definition 3.1), i.e. 

there exists a Borel subset V C U such that p(X \ V) = 0. (1) 

We are tempted to say that p is carried by U , but at this point it is not assumed that 
U is a Borel subset of X. This allows the framework to be applied to systems for which 
it is not known whether the subset of solutions is a Borel set in the appropriate space. 
Nevertheless, U is certainly measurable with respect to the Lebesgue completion p of 
p, and p(U) = 1 (see Remark 3.2). 

The terminology of trajectory statistical solution is inspired by the notion of trajec¬ 
tory attractor given, for instance, in [18] (see also [65]), and the connection between 
the two will become more apparent in a future work, when we focus on the case of a 
stationary trajectory statistical solution. 

For the initial value problem, we consider the interval I as being closed and bounded 
on the left, with left endpoint t 0 , and consider a tight Borel probability measure /io, 
defined on the space X, as the initial probability distribution of the state of the 
system. We also consider the projection operator Il io : X —> X which takes u , in X, 
into its value h [ to u = u(t 0 ), at time to- Then, the initial value problem for trajectory 
statistical solutions is simply to find a ^/-trajectory statistical solution p such that 
IR 0 p = /i 0 (see Problem 3.1). 

The existence of trajectory statistical solutions for this initial value problem relies 
on essentially three conditions: (i) a set-theoretic surjective condition on the pro¬ 
jection of U at the initial time, which in the applications is simply a statement of 
existence of “global” solutions, over the whole interval /, for every initial condition 
in the “phase” space X ; (ii) a topological compactness condition on the subset U 
associated with a certain family of compact sets of initial conditions, which in ap¬ 
plications follows from appropriate compact embedding theorems typically used for 
the existence of individual solutions of the system; and (iii) a joint topological and 
measure-theoretic property saying that the family of compact sets of initial condi¬ 
tions in the previous condition is sufficiently large to approximate from below any 
tight Borel probability measure on X (see Theorem 3.1). 

It might happen that the set of initial conditions for the existence of individual 
solutions does not coincide with the space in which the continuity in time holds. In 
this case, the result can in fact be modified to yield statistical solutions only carried 
by the “good” set of initial conditions (see Theorems 3.2 and 3.5 and Remark 3.7). 

Next, we turn to the concept of statistical solution in phase space. At this point, we 
need to be more specific about the evolution problem. Hence, we look for statistical 
solutions of an evolution equation of the form 

u t = F(t,u). 


( 2 ) 
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For this differential equation to make sense, we need a vector space structure. With 
that in mind, we start with the same setting as above and add another Hausdorff 
space Z and a topological vector space Y. such that Z C X C Y^„, with continuous 
injections, and where Y ((.* is the dual of Y endowed with the weak star topology. We 
denote the duality product between Y and Y' by (•, -)y\y- We consider a function 
F : / x Z —)• Y' and assume that u G U is such that u{t) belongs to Z for almost 
every t G /. With this framework, the definition of statistical solution in phase space 
is that of a family of measures {pt}tei which satisfies (2) in a suitable weak-star sense 
in the mean, i.e. 

J $(u) dp t (u) = J {F(t,u),&(u)) Y ',Y dn t (u), (3) 

in the distribution sense on /, for appropriate cylindrical test functions $ : Y' —> R, in 
conjunction with some measurability and integrability properties (see Definition 3.2). 
If a statistical solution in phase space {pt}tei is obtained as the family of projections 
p t = h l t p, t G /, of a trajectory statistical solution p, we say that it is a projected 
statistical solution. 

For the initial value problem for statistical solutions in the phase space, we first 
make the same assumptions as above for the subset U , hence obtaining a W-trajectory 
statistical solution starting with a given initial measure pq. Then, we assume that, 
for every u £ W, the function ( u(t),v)Y\Y is absolutely continuous on /, for every 
v G Y, with (2) being satisfied in a weak sense, namely that 

^{u(t),v) Y ',Y = {F(t,u(t)),v) Y ',Y, (4) 

in the distribution sense on /, for every v G Y. We add measurability conditions on the 
spaces and on F which guarantee the appropriate measurability of the Nemytskii-type 
operator (£, u) H > F(t, u(t )). One measurability condition is that every Borel subset of 
Z be also a Borel subset of X, which is a condition satisfied in all the examples that we 
are aware of (see Section 2.3). We also add integrability conditions on the right hand 
side of (4), for each trajectory in U , and an associated integrability condition on the 
initial measure. Those integrability conditions are related to a priori estimates which 
are natural for the system. Then, under these conditions, we obtain that the initial 
value problem has a statistical solution in phase space for any tight Borel probability 
measure satisfying the integrability condition for the initial measure (Theorem 3.4). 
Part of the proof of this result consists in showing that, under suitable conditions, the 
family {p t } te j of the projections p t = Il t p, t G I, of a trajectory statistical solution p 
is a statistical solution in phase space (see Theorem 3.3). 

The result on the existence of a statistical solution for the associated initial value 
problem is complemented with a result saying essentially that any energy-type in¬ 
equality (or equality) valid for the individual solutions holds also, in average, for the 
trajectory statistical solution and, hence, also for the projected statistical solutions 
(Propositions 3.1 and 3.2). 
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The proof of existence of a trajectory statistical solution for the initial value prob¬ 
lem is based on the Krein-Milmam approximation of the initial measure by convex 
combinations of Dirac deltas, as done in [33, 36]. We construct a net of trajectory 
statistical solutions with the measure at the initial time approximating the given ini¬ 
tial measure. At the limit, we obtain the desired trajectory statistical solution. In 
order to pass to the limit, one needs a compactness result for measures. A suitable 
result of this kind in our abstract framework is the one developed by Topsoe along 
his work on a generalization of Prohorov’s Theorem ([57]) to spaces which are not 
necessarily Polish ([71, 72, 73]). 

Topsoe’s topology is based on semi-continuity and is finer than the weak-star topol¬ 
ogy, hence it has less compact sets, but compactness is not the main problem here. 
Instead, the finer topology is needed to yield more open sets, so that two different 
tight measures can be separated by open sets. As a consequence, the space of tight 
measures on X under this semi-continuity-weak-star topology is a Hausdorff space, a 
fact that does not hold in general for the weak-star topology. If, however, the space 
X is also completely regular, then both topologies do coincide. (See Section 2.4.) 

After completing the abstract theory, we present some applications of our frame¬ 
work. The first and natural one is the system of Navier-Stokes equations, on which 
our whole abstract formulation was based (Section 4.1). The second one is a reaction- 
diffusion equation (Section 4.2). The third and last example is a nonlinear hyperbolic 
wave equation (Section 4.3). In each case, the system is formulated in the abstract 
framework and the hypotheses that yield both trajectory statistical solution and 
statistical solution in phase space for the corresponding initial value problems are 
verified. In all these examples, the spaces Z and Y are separable Banach spaces. In 
a future work, we plan to show how the theory applies to finding space-homogeneous 
statistical solutions for equations such as the Navier-Stokes equations on the whole 
space R 3 , in which Y will be taken as the inductive limit of a sequence of separable 
Frechet spaces. 


2. Basic Tools 

In this section, we introduce the basic concepts underlying our results. 

2.1. Function spaces. When working with measures on topological spaces, the topo¬ 
logical structure is of fundamental importance. In this regard, we recall a few concepts 
that play an important role in this work. Besides the fundamental notion of a Haus¬ 
dorff space, which is a topological space in which two distinct points can be separated 
by disjoint open sets, we recall that a topological space is completely regular when 
every nonempty closed set and every singleton disjoint from it can be separated by 
a continuous function. A completely regular space in which every singleton is closed 
is called a Tychonoff space. A topological space is said to be completely metriz- 
able when there exists a metric compatible with the topology of the space and under 
which the space is complete. A topological space is called Polish when it is separable 
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and completely metrizable. Separable Banach spaces and separable Frechet spaces 
are examples of topological vector spaces which are Polish spaces. 

When X is a topological vector space, we denote its dual by X' and the duality 
product is denoted by (•, -)x',x- When X is endowed with its weak topology, we denote 
the space by X w . Similarly, we consider X' endowed with the weak-star topology, in 
which case we denote it by X 4*. Notice that, for any topological vector space X, the 
space X' w * is always a Hausdorff locally convex topological vector space ([26, Section 
1.11.1]). If X is a Banach space, the norm in X is denoted by || • ||x, while || • ||x' 
denotes the usual operator norm in the dual space. 

Let X be a Hausdorff space and / clan arbitrary interval. Denote by C(I,X) 
the space of continuous paths in X defined on /, i.e. the space of all functions 
u : / —y X which are continuous. The compact-open topology in C(I,X) is the 
topology generated by the subbase consisting of sets of the form 

S(J, U) = {ue C(I, X) | u(J) C U}, 

where J is a compact subinterval of / and U is an open subset of X. When endowed 
with the compact-open topology, this space is denoted by X = Ci oc (J, A") and is a 
Hausdorff space. 

The subscript “loc” in C loc (7, X) refers to the fact that this topology considers com¬ 
pact sets in /. When X is a uniform space, the compact-open topology in C loc (J, A") 
coincides with the topology of uniform convergence on compact subsets ([45, 
Theorem 7.11]). This holds, in particular, when X is a topological vector space, which 
is the case in the applications that are presented in Section 4. 

For any t G /, let : X —* X be the “projection” map at time t defined by 

Yl t u = u(t), Vu G X. (5) 

It is readily verified that n f is continuous with respect to the compact-open topology. 

We also consider the space of bounded and continuous real-valued functions on X, 
denoted by C^(X). When X is a subset of IT 71 , m G IN, we also consider the space 
C£°(X) of infinitely differentiable real-valued functions on X which are compactly 
supported in the interior of X. 

2.2. Elements of measure theory. Let J be a topological space and %$x denote 
the cr-algebra of Borel sets in X. We denote by A4 (X) the set of finite and nonnegative 
Borel measures on X, i.e., the set of nonnegative measures /i defined on *3 X such that 
n(X) < oo. The subset of A4(X) consisting of Borel probability measures is denoted 
by V(X). The space Xi(X) can be identified with a subset of the dual space C^(X)' 
of the space C^(X). 

A carrier of a measure is any measurable subset of full measure, i.e., such that 
its complement has null measure. If C is a carrier for a measure /i, we say that /i is 
carried by C. If a probability measure is carried by a single point iGl, then it is a 
Dirac measure and it is denoted by 8 X . A probability measure that can be written 
as a (finite) convex combination of Dirac measures is called a discrete measure. 
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Given a family of sets $ C 23 x, we say that a Borel measure /i on X is inner 
regular with respect to the family ^ if 

n{A) = sup{/i(F) | F e £ and F C A }, \/A e 23x- (6) 

We say that a Borel measure /i is outer regular with respect to the family S if 

/i(A) = inf {fi(F) | F e S and A C F}, VA e 23 x . (7) 

A tight measure is a nonnegative Borel measure which is inner regular with respect 

to the family of compact subsets of X (such a measure is also called a Radon measure, 
see [9]). If a finite Borel measure n on A" is both tight and outer regular with respect 
to the family of open sets of X, then we say that /r is a regular measure. When X 
is a Polish space, every finite Borel measure is regular ([1, Theorem 12.7]). In case 
X is just a metrizable space, every finite Borel measure is inner regular with respect 
to the family of closed subsets of X and outer regular with respect to the family of 
open sets of X ([1, Theorem 12.5]) (such a measure is called normal in [1]). 

For a compact and metrizable space X , it follows in particular from the result in [1, 
Theorem 12.5] that every finite Borel measure is tight. The metrizability is indeed a 
necessary condition, since it is possible to construct a finite Borel measure defined on 
a certain no n m etrizable compact Hausdorff space which is not tight (see [1, Example 
12.9]). 

Furthermore, a net {p a }a of measures in Af(X) is said to be uniformly tight if 
for every e > 0 we can find a compact set K c X such that 

Ha(X\K) < e, Vo. 

The set of measures /i e A4(X) which are tight is denoted by At(X, tight). The 
subset of At (X, tight) consisting of probability measures is denoted by V(X, tight). 

Now consider a Hausdorff space Y and let F : X —y Y be a Borel measurable 
function. Then for every measure /i on 23 x we define a measure Ffi on 23 y by 

Ffi(E) = /i(F _1 (F)), VF e 23 r , 

which is called the induced measure from /i by F on 23y, also known as push- 
forward of /i by F. When fi is a tight measure and F is a continuous function, the 
induced measure F/i is also tight. 

In regard to the concept of induced measures, if (p : Y —y R is a F/i-integrable 
function then <p o F is /j-integrable and 

j tp o Fd/i = / ipdFfj, (8) 

J X JY 

(see [1, Theorem 13.46]). 

For the sake of notation, if /i e Af(X) and / is a /i-integrable function, we write 
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In the case of real numbers, we are also interested in the Lebesgue measure, which 
we denote by A, and in the Lebesgue subsets of intervals / C R. We denote the 
a-algebra of those sets by £/. 

2.3. Continuous injection of Borel sets. In the case of statistical solutions in 
phase space, in the abstract framework that we consider, a fundamental property 
that we need concerns the continuous injection of Borel subsets of a topological space 
into another topological space. More precisely, when we have two topological spaces 
Z and X, with Z continuously injected into X , meaning that there exists a continuous 
injective map j : Z —y X, we are interested in knowing whether the Borel subsets of 
Z are taken into Borel subsets of X by the injection j. Of course, this is equivalent 
to asking that the open subsets of Z are taken into Borel subsets of X. 

In general, this is a delicate issue. In fact, one can take X = [0,1] with the usual 
norm inherited from R and Z = [0,1] with the zero-one norm (associated with the 
discrete topology), so that every subset of Z is an open set, hence Borel in Z, but 
they are certainly not all included in the family of Borel subsets of X. It is worth 
noticing in this case that X is a compact metric space and Z is a locally compact 
metric space, but Z is not separable. 

In the case that Z and X are Polish spaces, then [9, Theorem 6.8.6] guarantees 
that j(B) is Borel in A", for any Borel B in Z. The theorem actually allows X to 
be more generally a Souslin space, which is defined as a continuous image of a Polish 
space. A closely related result can be deduced from [10, Lemmas 6.7.6 and 6.7.7], 
allowing X to be simply a metrizable space, while assuming that Z is a Lusin space, 
which is a space intermediate between Polish and Souslin, that can be characterized 
as the image of a Polish space under a continuous bijective map [10, Definition 6.4.6 
and Proposition 6.4.12]. Thus, if Z is a Lusin space and X is a metrizable space, 
then every Borel subset of Z is a Borel subset of X. This situation encompasses 
many important applications in which Z and X are Sobolev spaces, or other classical 
Banach spaces like Besov and Morrey, as long as Z is separable. 

Another particular situation important for us is when we consider the weak and 
strong topologies of a topological vector space. Since the strong topology is finer 
than the weak topology, every Borel set in the weak topology is also a Borel set in the 
strong topology. Conversely, if the topological vector space is separable and locally 
convex, then every strongly open set can be written as a countable union of strongly 
closed convex sets, and, thanks to the Hahn-Banach Theorem, every strongly closed 
convex set is weakly closed, so that every Borel set for the strong topology is also a 
Borel set for the weak topology. Therefore, in the case of separable locally convex 
topological vector space, both strong and weak Borel er-algebras in fact coincide. 
Combining this result with that of [10] mentioned above, we see that if Z is a Lusin 
space and X is obtained as a metrizable and separable locally convex topological 
vector space endowed with its weak topology, then the Borel subsets of Z are Borel 
subsets of X. This includes the case in which A is a separable Banach space, or even 
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separable Frechet, endowed with its weak topology, which is the case we consider in 
the applications in Sections 4.1 and 4.3. 

The results above can be easily extended to an important situation in partial differ¬ 
ential equations, namely when Z is the space of inhnitely-differentiable test functions 
with compact support and endowed with the topology used in the theory of distribu¬ 
tions, which is the inductive limit of a countable family of separable Frechet spaces 
[78, Section 1.1], Notice that any separable Frechet space is a Polish space. Now, if 
Z = U i^Zi is any inductive limit of a countable sequence of subsets Z % which are 
separable Frechet spaces, then any open set O in Z can be written as O = Uj(OflZj), 
and each O D Z, is open in Z,. Then, if X is either a metrizable space or a separable 
Frechet space endowed with its weak topology, we have that each Z, is also continu¬ 
ously injected into X, and, by the previous results, each O fl Z { is Borel in X. Thus, 
the countable union O = Uj(0 fl Zi) is also Borel in X. Hence, we deduce that any 
Borel subset of Z is a Borel subset of X. 

Finally, we should mention that the results described above can be extended to 
include products of those spaces, allowing us to tackle systems of equations. 

The conditions on Z and X described above certainly do not exhaust all the possible 
situations (see [62] for further conditions). For this reason, our main results on the 
existence of statistical solutions in the phase space (Theorems 3.3, 3.4, and 3.5, and 
Proposition 3.1) do not assume any extra structure on the spaces Z and A", instead 
just assume that Borel subsets of Z are Borel subsets of X, leaving the verification 
of this condition to the applications. 

2.4. Topologies for measure spaces and related results. In [72], Topsoe con¬ 
sidered a topology in A4( X) obtained as the smallest one for which the mappings 
/i ha p(/) are upper semicontinuous, for every bounded and upper semi-continuous 
real-valued function / on X. Topsoe calls this topology the “weak topology”, but in 
order to avoid any confusion we call it here the weak-star semi-continuity topol¬ 
ogy on A4(A). When a net {/x a } a converges to /i with respect to this topology, we 

.. wsc* 

write ji a —^ /i. 

A more common topology used in A4(A) is the weak-star topology, which is 
the smallest topology for which the mappings /i ha /i(/) are continuous, for every 
bounded and continuous real-valued function / on X. If a net {fi a }a converges to /i 
with respect to this topology, we denote fi a —^ /i. 

It is not difficult to see that the weak-star semi-continuity topology is equivalent to 
the topology obtained as the smallest one for which the mappings /i ha p(/) are lower 
semi-continuous, for every bounded and lower semi-continuous real-valued function 
/ on X. Then, considering now / as a bounded and continuous real-valued function 
on A, it follows that, for every a, b E R, the set 

{/i G M(X) | p(/) e (a, b )} = 

{p e M(X) | fi(f) e (a, Too)} n {/i e M(X) | n{f) e (-oo, 6)} (9) 
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is open in the weak-star semi-continuity topology, since / is in particular a lower and 
upper semi-continuous function. This shows that the weak-star topology is in general 
weaker than the weak-star semi-continuity topology. Moreover, according to Lemma 
2.1 below, if X is a completely regular Hausdorff space, then these two topologies 
coincide when restricted to the space M(X, tight). 

The following lemma summarizes the properties about the weak-star semi-continuity 
topology that we have mentioned and provides some additional useful characteriza¬ 
tions (see [72, Theorem 8.1]). 

Lemma 2.1. Let X be a Hausdorff space. For a net {p a } a i n M(X) and p e M.(X), 
consider the following statements: 

wsc* 

Fa ^ F> 

(2) limsup Fa(f) < F(f)> f or bounded upper semicontinuous function f; 

(3) liminf p a (f) > p(f), for all bounded lower semicontinuous function f; 

(4) lim a p a (X) = p(X) and limsup p a {F) < p(F), for all closed set F C X; 

(5) lim Q p Q (X) = p(X) and liminf p a (G) > p(G), for all open set G C X; 

w t 

(6) p Q ->■ p. 

Then the first five statements are equivalent and each of them implies the last one. 

Furthermore, if X is also completely regular and p e M(X, tight), then all six 
statements are equivalent. 

Although our framework is based on a general Hausdorff space, our proofs rely 
on reducing some structures to compact subsets, hence completely regular Hausdorff 
spaces, i.e., Tychonoff spaces. Moreover, since the measures in our proofs are usually 
tight, then in this setting both topologies coincide, so that we could have very well 
considered only the weak-star topology. However, we prefer to use the weak-star 
semi-continuity topology since it is a more natural topology for arbitrary Hausdorff 
spaces which simplifies our presentation and parts of the proofs, and which yields a 
compactness result in a stronger topology. 

When dealing with convergent nets in a given space, a natural question arises as 
to whether the limits are unique. This requires the given space of measures to be 
Hausdorff. The delicate issue is to determine the minimal hypotheses for that. 

If X is a metrizable topological space, then M.(X) is a Hausdorff space with respect 
to the weak-star topology (see [1, Section 15.1]), and hence also with respect to the 
weak-star semi-continuity topology. However, requiring X to be metrizable is too 
restrictive for our purposes. In looking for a more general setting for the space X, we 
were led to work within the space of tight measures A4(X, tight) and with the weak- 
star semi-continuity topology, which Topsoe proved to be a Hausdorff space. This 
key result was in fact a motivation for Topsoe to advance his work on the subject 
(see [72, Preface]). A proof of this fact is given in [72, Theorem 11.2] by showing that 
the limits of convergent nets are unique. Here we chose to include a proof showing 
directly that distinct measures in Ai(X) can be separated by open sets. 
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Theorem 2.1. Let X be a Hausdorff space. Then, Xi(X, tight ) is a Hausdorff space 
with respect to the weak-star semi-continuity topology. 


Proof. Consider two distinct measures pi, p 2 in Ai(X, tight). 

Suppose at first that /ii(df) ^ /r 2 (A"). Let us assume, without loss of generality, 
that lii(X) < /z 2 (X). Then there exists £ > 0 such that 

Ti(X) < Ti(X) + £ < /a 2 (X) — e < n 2 (X). (10) 

Denote a = gi(X) + £ and b = /j, 2 (X) — £. Let / be the constant function / = 1 on 
X. Then, it follows from (10) that 

/A e {/i G M(X) | n(f) G (-oo, a)} (11) 

and 

/r 2 g {/r g M(X) | g ( b , +oo)}. (12) 

Since / is in particular a bounded and continuous real-valued function on X, it follows 
that the sets in (11) and (12) are open sets in X. These sets are also clearly disjoint. 
Hence, pi and p 2 can be separated by disjoint open sets in X. 

Now suppose that ^(X) = p. 2 (X). We claim that there exists a set A e 23x such 
that /ii(H) < /j. 2 (A), where A and A denote the closure and interior of A, respectively. 
Indeed, suppose by contradiction that 

fH(A) >/! 2 (i), VH e 23 a -. (13) 


Consider E G 23x and let K \, K 2 be arbitrary compact sets in X satisfying K\ c 
X\E and K 2 C E. Then, since X is a Hausdorff space, there exist disjoint open sets 
/ j i . B 2 in X such that Ad C B \ and K 2 C B 2 . 

In particular, using that B 2 is an open set, it follows from (13) that 

ia 2 (B 2 ) = p 2 (B 2 ) < /ii(H 2 ). (14) 

But clearly B 2 C X\K\. Then, 

/ii(H 2 ) < /j,i(X\Ki) = mix) - m(Ki). (15) 

From (14) and (15), we obtain that 

H2{k 2 ) < /i 2 (H 2 ) < mix) - m(Ki). 


Now since K\ and K 2 were chosen arbitrarily, taking the supremum over all compact 
sets Ad, K 2 with Ad C X\E and K 2 C E, and using that (i\ and /j 2 are tight, it follows 
that 


/i 2 (£) < Ti(x) - m(X\E) = m(E). 


Thus, 


vag 2S x . 


(16) 


This implies in particular that 


m(x) - m{E) = p2{x\e) < m{x\E) = m{x) - bi(e). 
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Using the hypothesis that /ii(A") = p 2 ( X ), we then obtain 

fi 2 (E) > pi (E), VE e (17) 

Now (16) and (17) yield pi = p 2 , which is a contradiction. 

Thus, we may consider A e 23X such that pi(A) < /i 2 (H). The argument now 
follows analogously to the previous case. Consider e > 0 satisfying 

Pi(^4) < Pi(^4) + £ < P 2 (^4) — £ < P2(^4) 

and denote a' = pi(A) + £ and b' = p 2 (A) — £. Since the characteristic function xa 
of A is bounded and upper semi-continuous and the characteristic function xa °f 71 
is bounded and lower semicontinuous, then the sets 

{p e M(X) | h(xa) e (—oo,a')} 

and 

{p 6 M(X) | p(xa) e (&', +oo)} 

are clearly disjoint open sets in X containing p x and p 2 , respectively. Thus, pi and 
p 2 can also be separated by disjoint open sets of X in case fii(X) = p 2 (A"). 

This proves that M. (X, tight) is a Hausdorff space. □ 

Moreover, if X is assumed to be a completely regular Hausdorff space, then by 
Lemma 2.1 the weak-star semi-continuity and weak-star topologies are the same in 
M(X, tight). Thus, it follows from Theorem 2.1 that A4(X, tight) is also Hausdorff 
with respect to the weak-star topology. We have just proved the following corollary: 

Corollary 2.1. Let X be a completely regular Hausdorff space. Then A4(X, tight) is 
a Hausdorff space with respect to the weak-star topology. 

Remark 2.1. If we relaxed the hypothesis in Corollary 2.1 by assuming X to be a 
regular space instead of a completely regular space, then this result would no longer 
be valid. Indeed, it is possible to construct examples of regular spaces which are not 
completely regular and containing two distinct points a, b for which every continuous 
real-valued function satisfies /(a) = f(b) (see [55, 59]). The corresponding space 
M(X, tight) is then not Hausdorff with respect to the weak-star topology, since it 
suffices to consider the Dirac measures S a and S & concentrated on a and b, respectively, 
and to note that 5 a (f) = dfff), for every bounded and continuous real-valued function 
/ on X. 

We next state a result of compactness on the space of tight measures M(X, tight) 
that is essential for our main result. For a proof of this fact, see [72, Theorem 9.1]. 

Theorem 2.2. Let X be a Hausdorff space and let {/i a } a be a net in A4(X, tight) 
such that limsupp a (X) < oo. If {p a }a is uniformly tight, then it is compact with 
respect to the weak-star semi-continuity topology in M(X, tight). 
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The previous theorem allows us to obtain a convergent subnet of a given net in 
M(X, tight), provided it satisfies the required conditions. Also, by using Theorem 
2.1, we can guarantee that the limit of this convergent subnet is unique. 

Evidently, all the results shown above are also valid in the space of probability 
measures. In what follows, we consider both the spaces of probability measures 
defined over the Hausdorff space X and over the space of continuous paths X. 

2.5. Weak-star scalarwise derivative. For an evolutionary differential equation, 
we need a notion of time-derivative for functions with values in a vector space. For 
our purposes we only need a scalarwise-type derivative, but since we look for solutions 
which are also weakly continuous (see Section 2.6), this turns out to be a functional 
derivative associated with the functional integral in the sense of Pettis (see Remark 


2 . 2 ). 


We follow the nomenclature of scalarwise property given in [11, 26], except that we 
restrict ourselves to functionals with values in the dual of a topological vector space. 
More precisely, we consider a topological vector space Y and its dual Y', and we say 
that u : / —y Y’ has a certain property P weak-star scalarwise if each function 
t i — ^ ( u(t),v)Y',Y has the property P, for every v G Y. We will be interested in the 
weak-star scalarwise versions of the properties of a function being (Lebesgue) mea¬ 
surable, integrable, locally integrable, continuous, absolutely continuous, and almost 
everywhere differentiable. 

Saying that u is weak-star scalarwise continuous is the same as saying that u is 
weak-star continuous (i.e u is continuous as a function from / into the topological 
vector space Y^*). But the same cannot be said about the other properties. For 
instance, u can be weak-star scalarwise measurable without being Borel measurable 
from / into Y^, (see Remark 2.5). 

We are interested in functions u : I —> Y' that are weak-star scalarwise continuous 
and such that there exists a weak-star scalarwise locally integrable function w : I —> 
Y' which is the weak-star scalarwise derivative of u in the distribution sense on /, i.e. 



This is equivalent to assuming that u : I —> Y' is weak-star scalarwise continuous 
and that there exists a weak-star scalarwise locally integrable function w : I —> Y' 
satisfying 



This is also equivalent (see e.g. [37, Theorem 3.35]) to considering functions that 
are weak-star scalarwise absolutely continuous and such that there exists a weak-star 
scalarwise locally integrable function w : I —>■ Y' satisfying 



Vu G Y, for a.e. t G I. 


( 20 ) 
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Using the shorter property (20) of the equivalent forms (18), (19), and (20), we 
define the space 

u is weak-star scalarwise absolutely continuous and 
{{ . j ^ y; . 3 w : / -I h' weak-star scalarwise locally integrable 

with — {u(t), v)y>,y — ( w (t), v )Y',Y, a.e. t E /, Vn G Y 

CJ .L 

The weak-star scalarwise derivative w of an element of 3d, as given in the dehnition 
(21), is uniquely defined in the case that Y w is separable (see Remark 2.3), but may 
not be unique in general (see Remark 2.4). 

Remark 2.2. If w : I —> Y' is weak-star scalarwise integrable on /, then w induces a 
linear map A w belonging to the algebraic dual of Y, associating each element v G Y 
to the scalar A w (v) = J I {w(t),v)yi > Y d£. This map is uniquely defined. If this 
map belongs to Y', then it is the Pettis integral of w over / and it is written as 
A w = fjVj(t) dt. Similarly for the locally weak-star scalarwise integrable case, when 
integrating over compact measurable sets within /. More generally, one may consider 
integrals on any topological vector space, not necessarily dual spaces, but this is the 
context we are interested in. There are a number of results concerning the scalarwise 
measurability and integrability and conditions for the existence of the Pettis integral; 
see e.g. [26, 11] for further details. For instance, if Y is a locally convex topological 
vector space with the so-called GDF property (i.e. any linear map of Y into a Banach 
space B having a graph that is sequentially closed in Y x B is a continuous map), then 
any (weak-star) scalarwise integrable function w from / into Y' has a Pettis integral 
well defined and belonging to Y' [26, Theorem 8.16.1]. Any Frechet space and any 
inductive limit of Frechet spaces have the GDF property [26, Subsection 8.14.12], In 
our work, however, we do not need to worry about conditions for the existence of 
the Pettis integral in Y 1 , since the weak formulation is all we need. Nevertheless, the 
Pettis integral of any weak-star scalarwise derivative of a function in 3d does exist 
since the equivalent formulation (19) shows that the weak-star scalarwise integral is 
a continuous linear function on Y. 

Remark 2.3. In the case that Y w is separable, i.e. if Y is separable under the 
weak topology, then w : I —> Y' satisfying the conditions in the dehnition of 3d 
in (21) is unique (up to measure zero on I). Indeed, suppose W\ and w 2 are two 
functions satisfying the conditions on w in the dehnition of 3d- Assume there exists 
a countable set D = { Vj}j in Y which is dense in Y with respect to the weak 
topology. Then, for each Vj, we have (uq (t),Vj)Y',Y — ( w 2 if)i v j)Y',Yi f° r al mos t every 
t in I. Since the set D is countable, there is a set of full measure Iq in / for which 
( w 2 (t) — Wi(t),Vj)Y',Y = 0, V£ G Jo, Vj G IN. For each t G Jo, since D is dense in Y in 
the weak topology, we find that (w 2 (t) — W\(t),v)y',Y = 0, for all v EY. This means 
that w 2 (t) = W\(t), for all t E Iq. Since Iq is of full measure, we find that W\ = w 2 up 
to a set of measure zero. 
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Remark 2.4. If Y w is not separable, then we cannot guarantee that w : I Y' 
satisfying the conditions in the definition of 34 in (21) is unique. Indeed, take any 
nonseparable Hilbert space Y with a basis with the cardinality of the continuum. 
Denote the inner product by (•> ')y and identify the dual Y' with Y itself. Let / = [0, 1 ] 
and let {u s } se / be an orthonormal basis, parametrized by the continuum I. Consider 
the function w : / —> Y given by w(t) = Vt, for t G [0,1]. For each v s , the real-valued 
function 1 H > ( w(t),v s )Y is equal to 0, for t ^ s, and to 1, for t = s. Thus, it is almost 
everywhere equal to zero and, in particular, Lebesgue integrable on I. If v G Y is 
arbitrary, then v is a linear combination of at most a countable set of indices, so that 
(■ w(t),v)Y is different from zero at most at a countable set of values t. This means 
that the real-valued function t i—>■ ( w(t),v)y is, again, zero almost everywhere and 
Lebesgue integrable. Thus, 

i 

(w(t),v)Y dt = 0, Vu G Y. 

Nevertheless, w is not zero almost everywhere in Y. In fact, it is different from zero 
everywhere on I. Thus, if a function w : I Y' satisfies the conditions in the 
definition of 34 in (21), then so does w + w, which shows the nonuniqueness of the 
weak-star scalarwise derivative in this case. 

Remark 2.5. The same function w defined in Remark 2.4 when Y is a nonseparable 
Hilbert space with an orthonormal basis with the cardinality of the continuum is an 
example of a function which is (weak-star) scalarwise Borel measurable on / but not 
Borel measurable from / into Y w (which in this case is exactly T/J. Indeed, for each 
f G I, let O t be a weakly open set in Y which only contains the vector v t of the basis, 
not v s , for s ^ t (say O t = {v E Y : (v,Vt)y > 1/2}). Take any nonmeasurable 
subset A of I. Then, O = U t ^AOt is a weakly open set in Y such that iu _1 (0) = A is 
not measurable. 

Remark 2.6. In case Y is a Banach space, then so is Y'. and one may consider 
the Bochner integral in Y' [78]. The set of locally Bochner integrable functions that 
have a weak derivative which is also locally Bochner integrable in Y' forms the space 
M4oc (/; H'), and, in this case, the time derivative u t of a function u G H 7 /,’/ (I]Y r ) 
is uniquely defined and satisfies u(ti) = u(t 0 ) + u t (t) dt in Y', for almost every 
to, ti G /, where the time integral is the Bochner integral in Y'. In this case, we notice 
that H 7 ]// (/; Y') is included in the space 34 defined in (21) (see e.g. [69, Lemma 3.1.1]), 
with (18), (19), and (20) holding with w = u t . Therefore, the framework above can 
be used in applications where the differential equation holds in a functional sense 
associated with the Bochner integral. 

2.6. Evolution equations. The statistical solutions in a phase space X (see Defi¬ 
nition 3.2) are directly related to an evolution equation of the form 

«t(<) = F(t,u(t)), 



( 22 ) 


ABSTRACT FRAMEWORK FOR THE THEORY OF STATISTICAL SOLUTIONS 


17 


where the unknown u belongs to the space X = C\ oc (I, X), with X being a Hausdorff 
space, and I C R, an interval. 

We want to be able to make sense of the equation (22) under minimal hypotheses 
on the structure of F and on the spaces involved. Of course, the peculiarities of the 
solutions are of interest in each application, but for our general result of existence 
of statistical solutions for the initial value problem, we want to avoid superfluous 
conditions. 

For partial differential equations, the right hand side of the equation (22) involves 
spatial derivatives of the solution u, so it is not expected that F be defined on X, but, 
instead, on a more regular space included in X. Moreover, for the time-derivative to 
be defined, we need a vector space structure, so that u needs to be included in some 
topological vector space, possibly less regular than X. 

With that in mind, we consider another Hausdorff space Z and a topological vector 
space Y, and consider F such that 

F:/xZaF, (23) 

with the spaces satisfying 

2clcr;„ (24) 

where each space is continuously included in the next one, and where Y' v , is the dual 
space Y' endowed with the weak-star topology. The reason for the inclusion X C Y^* 
is that we will in fact consider the weak formulation of the evolution equation (22), 
so that (24) is the natural assumption to make. In many cases, Z — Y , but this is 
not always convenient or necessary (see Section 4.2 and Remark 3.4). 

Then, we look at the weak formulation of (22), i.e. 

^{u(t),v) Y ', Y = {F(t,u(t)),v) Y ’,Y, (25) 

in the distribution sense on /, for any v G Y. 

Since u G X = Ci oc (/,X) and X C Y^», we have that u G C\ oc (I,Y ((.*), so that 
t i — x ('u(t), v)y',y is continuous, for every v G Y. Thus, in order for (25) to make 
sense, t H > (u(t),v)Y',Y needs to be absolutely continuous and the right hand side of 
(25) needs to be integrable. With that in mind, we consider u belonging also to the 
set 3d defined in (21). Moreover, we need u such that u(t) G Z for almost every tel. 
Therefore, we define the spaces 

Z = {u G X : u{t) G Z for almost all tel}, (26) 

and 

x 1 = zny 1 . (27) 

Then, if u e X i is such that 

1 1 —y F(t,u(t )) is weak-star scalarwise locally integrable, (28) 

as a function from / into Y', it makes sense to require that Ut = F(t,u ) be valid in 
the weak sense (25). Notice that, for u G X\ such that (28) holds, condition (25) is 
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precisely (18) with w(t) = F(t,u(t )) and is equivalent to both (19) and (20), with 
the same w. 

2.7. Cylindrical test functions. Consider a topological vector space Y and let 
vi,...,Vk G Y, where k G IN. Let 0 be a continuously differentiable real-valued 
function on R k with compact support. For each u G Y', define $(«) G R by 

$(u) = (j)({u, Vi) Y ',Y, •••,(«, V k ) Y ',Y)- 

The function is clearly continuous from Y' to R and in fact it is differentiable in 
Y', with differential at m G Y' given by 

k 

$'(«) = X v i)y',y, ■•■,(«, v k ) Y ’,Y)vj, (29) 

3 = 1 

where dj(f> denotes the derivative of <f> with respect to its j-tli coordinate. For each 
u, w G Y', 

k 

(w, &(u)) Y ' X = ^dj(j)((u,vi) Y \ Y , • • •, {u,v k ) Y ' iY ){wiVj) Y ' iY 

3 = 1 

is the Gateaux derivative of $ at u , in the direction of w. Moreover, $'(«) belongs to 
Y and the differential <f>' is continuous from Y', into itself, where we recall that Y' 
denotes the space Y' endowed with its weak-star topology. If Y is a Banach space, 
then Y' is a Banach space under the strong operator norm, and d*' is the Frechet 
differential of <f> in this strong norm of Y'. 

Functions of this form are called cylindrical test functions in Y' and play an 
important role as test functions in the definition of statistical solution in phase space 
(see Definition 3.2). In that context, we will also consider a Hausdorff space X which 
is assumed to be continuously imbedded in Y^. Notice that since u i —> (u,v) Y ' jY is 
continuous in the weak-star topology for any v G Y, the function $ is also continuous 
from Y' wif into R. Then, since X is continuously imbedded in Y' i;ir , we may consider $ 
restricted to X , which is continuous as a function from X into R. 

Remark 2.7. The set of cylindrical test functions is a relatively large set, as can 
be seen from the Stone-Weierstrass Theorem (see e.g. [25, Theorem IV.6.16]). In 
fact, consider a compact set K' in Y^, (e.g. a (strongly) closed ball in Y 1 , if Y is a 
normed space or, more generally, any closed subset of Y^ t which is weakly bounded 
and equicontinuous ([26, Theorem 1.11.4])). Denote by Sk' the collection of the real¬ 
valued functions on K' which are the restriction to K ' of the cylindrical test functions. 
Clearly, Sk' C C(K' ) and, if Ti, T 2 G Sk ', then their sum Ti + T 2 and their product 
T 1 'L 2 also belong to Sk'- This means that Sk' is a subalgebra of the algebra C(K'). 
Choose any v EY and notice that the set J = {(u 1 v) Y ^ Y : u G K'} is compact in R, 
otherwise K’ would not be weak-star compact. Then, by taking any 0 : R —)■ R which 
is continuously differentiable and is equal to 1 on J, we see that <F(n) = cp((u,v) Y ^ Y ) 
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is a cylindrical test function which is equal to 1 for u G K', showing that Sk' contains 
the unit element. Moreover, if Ui,U 2 are distinct points in K there exists v G Y 
such that (ui,v)y',y ,v)y',y- Hence, by choosing a continuously differentiable 

function 0 : R —* R which is compactly supported and assumes different values at 
the points ( u\,v)y>,y and (u2,v)y',y, we see that <h(u) = v)y>,y) is a cylindrical 

test function which assumes different values at u± and U 2 , proving that Sk> separates 
the points in K'. Therefore, the Stone-Weierstrass Theorem yields that Sk' is dense 
in C(K'), for the uniform topology. Similarly, since X is continuously imbedded into 
Y^, it follows that, for any compact subset K of X , the collection of the real-valued 
functions on K which are the restrictions to K of cylindrical test functions is dense 
in C{K). 

2.8. The Nemytskii operator. In relation to the evolution equation (22), we con¬ 
sider a function F of the form F : / x Z —y W, where Z and W are Hausdorff topolog¬ 
ical spaces, and consider the operator G : I x Z -A- W given by G(t,u ) = F(t,u(t)), 
where Z is defined in (26). It is of fundamental importance to deduce measurability 
properties for G from those of F. This is the aim of this section. This extended 
operator G is known as a Nemytskii operator in the context of partial differential 
equations, where the Hausdorff spaces above are simply spaces of real-valued func¬ 
tions defined on subdomains of Euclidean spaces. In this regard, we have the following 
results. 


Lemma 2.2. Let X and W be Hausdorff spaces and let I be an interval in R. Con¬ 
sider a function F : I x X —>■ W and define the associated Nemytskii operator 
G : I x X -A- W by G{t,u) = F(t,u{t)), for t e I and u e X = C loc (I,X). If 
F is a (£/ <E> Q3x, ^w)-measurable function, then G is (£/ ® Q3*, 23 w)-measurable. 


Proof. Consider the projection 


n 7 : I x X —>■ I 
(■ t,u ) i —y t, 

and the evaluation operator 

U : I x X -A X 

(t,u) i-A- u(t). (30) 

The projection n 7 is clearly a (£/ <g> 23^, £/)-measurable function. Moreover, U is 
a continuous function and then, in particular, (23/ x ;r, 53 x) -measurable. Since X 
is a Hausdorff space and / is a second countable Hausdorff space it follows that 
53/x.r = 53/ <E> 53^ ([9, Lemma 6.4.2]). Thus, U is also a (£/ (g) IB x, 53x)-measurable 
function. From this we obtain that the function (H/, U) : / x X -A- / x X , defined 
by (n j,U)(t,u) = (t,u(t)), is (£/ <g) 53 x,£i 53x)-measurable. Now, since F is, by 
assumption, (£/ (g) 53x, 23 w )-measurable, and G can be written as G = F o (n/, U ) 
we conclude that G is (£/ (g) fB^, 23 w) - measurable. □ 
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Now we consider the measurability of the trivial extension of a function F defined 
on / x Z to / x X. 


Lemma 2.3. Let Z, X, W be (nonempty) Hausdorff spaces such that Z C X, with 
continuous inclusion, and let / be an interval in It. Consider a function F : I x Z -E 
W. Let w be an arbitrary element in W and define the extension F : I x X —>■ W by 


F{t,u) 


F(t, u ), if u E Z, 
w, otherwise. 


(31) 


Assume that every Borel subset of Z is a Borel subset of X and that F is a (£/ <g) 
?Bw)-measurable function. Then F is a (£/ <E> Q3x, 53iy )-measurable function. 


Proof. Let E E and note that 

F-'(E)U(Ix(X\Z)), if weE, 

' \ ifw^E. 

Since F is (£/ <g) 53z, %$w) -measurable then F" 1 (L’) E £/ <8> From the hypothesis 
that *Bz C this implies that F" 1 (i?) E £/ <Z ^x- Moreover, again from *3 Z C 
53x, we have, in particular, that Z is Borel in A", so that X \ Z E 53x, as well, and, 
hence, I x (A \ Z) E £/ <g) Thus, F~ l (E) E £/ ® 5Sx, regardless of whether w 
belongs to E or not. This proves that F is (£/ ® 53x, 53vi/)-nieasurable. □ 


Combining Lemmas 2.2 and 2.3 we obtain the following result. 


Proposition 2.1. Let Z, X, W be (nonempty) Hausdorff spaces such that Z C X, 
with continuous inclusion, and let I be an interval in R. Consider a function F : 
I x Z W. Let w be an arbitrary element ofW and define the associated Nemytskii 
operator G : I x X —>■ W by 


G(t , u ) 


F(t, u(t)), if u(t) E Z, 
w, otherwise , 


(32) 


where X = Ci oc (J, A). Assume that every Borel subset of Z is a Borel subset of X 
and that F is a (£/ ® ^&z,^w)- measura ble function. Then G is (£/ <g) 53^,53vr)- 
measurable. 


3. Abstract Results 

In this section, we present our abstract framework for the theory of statistical 
solutions. First, in Section 3.1, we give our general definitions of statistical solutions 
in trajectory space and in phase space. Then, in Sections 3.2 and 3.3 we prove the 
main results on the existence of these general types of statistical solutions with respect 
to a given initial data. 
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3.1. Types of Statistical Solutions. We first define statistical solutions in the 
space of continuous paths X = C\ oc (I,X), in a Hausdorff space X. They are named 
trajectory statistical solutions, owing to the fact that they are measures carried 
by a measurable subset of a certain set U in X which, in applications, would consist 
in the set of trajectories, i.e. the set of solutions, in an appropriate sense, of a given 
evolution equation. At this abstract level, however, there is no evolution equation, 
and the problem is simply to find a measure carried by a given subset U of X. As 
such, this is a trivial problem, as showed in Remark 3.5. The interesting and difficult 
problem is the corresponding Initial Value Problem 3.1. Nevertheless, we start with 
the following definition. 

Definition 3.1. Let X be a Hausdorff space and let I C R be an arbitrary interval. 
Consider X = C\ oc (I,X) and let U be a subset of X . We say that a Borel probabil¬ 
ity measure p on X is a U-trajectory statistical solution over I (or simply a 
trajectory statistical solution) if 

(i) p is tight; 

(ii) p is carried by a Borel subset of X included in U, i.e., there exists V G 23^ 
such that VcW and p(X \ V) =0. 

Remark 3.1. Our abstract definition of a trajectory statistical solution was inspired 
by the concept of a Vishik-Fursikov measure given in [36]. Such measures are defined 
within the context of the Navier-Stokes equations and have the property of being 
carried by their set of weak solutions, called Leray-Hopf weak solutions. In [36, 
Propositions 2.9 and 2.12] it is proved that the set of Leray-Hopf weak solutions is a 
Borel set in the corresponding space of continuous paths. However, since we do not 
know whether this is always the case in every application, we prefer not to assume 
that U is Borel, and assume instead that there exists a Borel subset of U that carries 
the measure p. 

Remark 3.2. From the Definition 3.1, however, we see that U\V belongs to the 
Borel null set X \ V, hence we can certainly say that U is measurable with respect to 
the Lebesgue completion of p , which we denote by p, and so that p is carried by U. 

Now we define statistical solutions in phase space. Unlike the definition of a trajec¬ 
tory statistical solution, which is given by a single measure defined on X , this second 
type of statistical solution consists in a family of measures defined on the Hausdorff 
space X and parametrized by an index t varying in an interval / C R. For the sta¬ 
tistical solutions in phase space, we do need an evolution equation, although minimal 
hypotheses will be needed, as described in Section 2.6. We then have the following 
definition of statistical solution in phase space: 

Definition 3.2. Let X and Z be Hausdorff spaces and Y be a topological vector space 
such that Z C X C Yff, with continuous injections, where Y' denotes the dual space 
of Y and Yf» is the dual space Y' endowed with the weak-star topology. Let I C R be 
an arbitrary interval and consider a function F : I x Z Y'. We say that a family 
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{pt}t£i of Borel probability measures in X is a statistical solution in phase space 
(or simply a statistical solution) of the evolution equation u t = F(t,u), over the 
interval I, if the following conditions are satisfied: 

(i) The function 

t 1 t I ip(u) d p t (u) 

J x 

is continuous on I, for every p E C^(X). 

(ii) For almost every t E I, the measure p t is carried by Z and the function 
u i—^ ( F(t,u),v)Y',Y is pt-integrable, for every v eY. Moreover, the map 


1 1 —y 


(F(t,u),v) Y ',Y d Pt(u) 


belongs to L\ oc (I), for every v E Y. 

(iii) For any cylindrical test function $ in Y', it follows that 

[ $(u)dp t (u)= [ $0) dpt>(u) + [ [ (F(s,u),&(u)) Y ',y dp a (u)ds, 
J x J x 

for all t,t' E I. 


(33) 


/,' JX 


In Definition 3.2, equation (33) represents a Liouville-type equation similar to that 
from statistical mechanics. In order to motivate the definition (33), let us suppose that 
a particular statistical solution {p t }tei is given in the form of a convex combination 
of Dirac measures, 

1 N 

pt = teI > 
n =1 

with equal probability 1/N, where N E N, and each u n is a smooth solution (say 
continuously differentiable as a function from the time interval / into the space Y') 
of the system 

u t (t) = F(t,u(t)), tEl. 

We then formally have 


d 

ds 


$ («) dp*(u) = $ K( S )) 

n= 1 ' 


Y'Y 



n =1 


jf m ^( s )), ^(^( s )))^^ 

n =1 


= / (F(s,u),&{u))y',y dp s (u). (34) 


Jx 

Thus, integrating with respect to s on [t',t\ yields (33). 

In Section 3.3 (Theorem 3.3), we prove that the family of measures obtained as 
the projections of a trajectory statistical solution at each t E I on the space X is a 
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statistical solution. This shows that every trajectory statistical solution on X yields 
a statistical solution on X. However, the converse is not necessarily true. We then 
call a statistical solution for which the converse is valid, i.e., which can be written 
as the projections on X of a trajectory statistical solution, a projected statistical 
solution, as given below. 

Definition 3.3. Let X and Z be Hausdorff spaces and Y be a topological vector 
space such that Z C X C Yf,, with continuous injections, where Y' denotes the dual 
space of Y and Yf, is the dual space Y' endowed with the weak-star topology. Let 
I C R be an arbitrary interval and let U be a subset of X. Consider a function 
F : / x Z —)• Y'. We say that a family {pt}tei of Borel probability measures on X 
is a statistical solution projected from a U -trajectory statistical solution 
of the evolution equation u t = F(t,u), over the interval L (or simply a projected 
statistical solution), when {pt}tei is a statistical solution of the evolution equation 
u t = F(t,u ) in the sense of Definition 3.2 and there exists a U-trajectory statistical 
solution p such that p t = Ll t p, for every t G /. 

Remark 3.3. In Definition 3.2, condition (iii), since $'(w) is a linear combination of 
vectors in Y (see (29)), the condition (ii) is sufficient to guarantee that the second 
integral on the right hand side of (33) is well-defined. More details in the proof of 
Theorem 3.3. 

Remark 3.4. In most applications, the spaces Z and Y in the Definition 3.2 of 
phase-space statistical solution are taken to be the same. Nevertheless, this is not a 
requirement in the abstract setting and we allow for this flexibility. This situation is 
illustrated in Section 4.2. 

Remark 3.5. Note that whenever U is a nonempty set, we can always obtain a 
trajectory statistical solution by considering the Dirac measure S u . for any element 
u eW (6 U is tight and {n} is a Borel set in U satisfying h u ({n}) = 1). A statistical 
solution can then also be easily obtained by considering the family of projections 
{d u (t)}tei- However, our main concern is not simply the existence of a measure or 
a family of measures satisfying the properties described in Definitions 3.1 or 3.2, 
respectively. Our aim is to prove the existence of such solutions for an initial value 
problem, as described below. 

In the case of trajectory statistical solutions, the initial value problem takes the 
following form: 

Problem 3.1 (Initial Value Problem for Trajectory Statistical Solutions). Let / C R 
be an interval closed and bounded on the left, with left end point t 0 , and let X be a 
Hausdorff space. Let X = C\ oc (L,X) be the space of continuous paths in X endowed 
with the compact-open topology. Let U be a given subset of X. Given an “initial” 
tight Borel probability measure po on X, we look for aL/-trajectory statistical solution 
p on X satisfying the initial condition n io p = p 0 , i.e., we look for a measure p G V(X) 
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satisfying conditions ( i ) and (ii) of Definition 3.1 and such that 

/) (n t - 0 1 0)) = MoO), 

The corresponding problem for statistical solutions is stated analogously: 

Problem 3.2 (Initial Value Problem for Statistical Solutions). Let X and Z be 
Hausdorff spaces and Y be a topological vector space such that Z C X C Yf t , with 
continuous injections. Let / C R be an interval closed and bounded on the left, 
with left end point t 0 . Consider F : I x Z Y'. Given an “initial” tight Borel 

probability measure p 0 on X , we look for a statistical solution {pt}t&i of the equation 
u t = F(t,u), over the interval /, and satisfying the initial condition p to = po- 

Note that in the Initial Value Problems 3.1 and 3.2 the interval / is considered as 
being closed and bounded on the left, with left end point t 0 . This point to represents, 
of course, the initial time in an application, so that po is the initial measure. 

3.2. Existence of Trajectory Statistical Solutions. In order to obtain the exis¬ 
tence of trajectory statistical solutions in the sense of Definition 3.1 and satisfying a 
given initial data (Problem 3.1), we impose some natural assumptions on the subset 
U C X. In fact, our main existence result for the initial value problem for trajectory 
statistical solutions reads as follows. 

Theorem 3.1. Let X be a Hausdorff space and let I be a real interval closed and 
bounded on the left with left end point t 0 . Let U C X be a subset of X = Ci oc (J, X). 
Suppose that 

(Hi) n to u = x. 

Suppose, moreover, that there exists a family .h'(A) of compact subsets of X such that 
(H2) Every tight Borel probability measure po on X is inner regular with respect to 
the family .h'(X) in the sense of (6); 

(H3) For every K e $!(X), the subset Llf'K D U is compact in X. 

Then, for any tight Borel probability measure po on X there exists a U-trajectory 
statistical solution p on I such that n to p = p 0 . 

Before going into the proof of Theorem 3.1, let us first discuss the hypotheses of 
the theorem and the main ideas in its proof. 

In the applications to a given evolution equation, the Hausdorff space X plays 
the role of the phase space associated with the equation, and the set U is the set of 
solutions in a given sense. These solutions are assumed to be continuous functions 
defined on a real time-interval / and with values in the phase space X , so that U is 
a subset of X = C\ oc (I,X). 

Then, hypothesis (HI) is simply a statement about global existence of solutions 
for every initial condition in A". Indeed, it requires that for every initial condition 
Mo 6 X, there is a solution u 6W, existing for the whole time interval /, and satisfying 
u(t 0 ) = u 0 . 
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Hypotheses (H2) and (H3) come together and are essentially a compactness condi¬ 
tion on the subsets of solutions with initial conditions in sets belonging to a certain 
family of compact subsets of the phase space which is sufficiently large to approxi¬ 
mate from below every tight Borel probability measure on A". In many applications, 
.ft'(A) may be considered as the entire family of compact sets of A, so that hypothesis 
(H2) holds trivially, and then hypothesis (H3) follows from usual a priori estimates 
and some compactness embedding theorem. This is the case, for instance, for the 
reaction-diffusion equation presented in Section 4.2. 

In other applications, however, such as to the incompressible Navier-Stokes equa¬ 
tions and the nonlinear wave equation (Sections 4.1 and 4.3), the compactness of 
n^K nu is not known to hold for every compact set K in A. The reason is related 
to the fact that the solutions are taken to be weakly continuous in time while being 
strongly continuous at the initial time (with the phase space A being equal to a sep¬ 
arable Banach space, in a set-theoretic sense, but endowed with the corresponding 
weak topology). It is precisely this strong continuity at the initial time which might 
be lost for the limit points of n^ 1 /! D U. In order to overcome this problem, the 
family IT (A) is taken to be slightly smaller, made only of the strongly compact sets. 
Using the energy inequality, one proves that the strong compactness of K implies that 
the limit points of nD U are also strongly continuous at the initial time, so that 
Tl^K ft IA is, in fact, compact in X . On the other hand, since in a separable Banach 
space the Borel cr-algebras corresponding to the strong and the weak topologies coin¬ 
cide (see Section 2.3), any Borel probability measure on A (endowed with the weak 
topology) is also a Borel measure with respect to the strong topology. Moreover, since 
every separable Banach space is also a Polish space, any Borel probability measure 
with respect to the strong topology of X is inner regular with respect to the family 
of (strongly) compact subsets (see Section 2.2), showing that hypothesis (H2) also 
holds. 

Let us now outline the main ideas of the proof itself. Starting with an initial 
measure /xo in V(X, tight), at a given time to, our intention is to show the existence 
of a measure p which is a trajectory statistical solution satisfying the initial condition 
n t 0 p = po- As usual, this measure p is obtained from the limit of a convergent net of 
measures. 

We first consider the case when the initial measure /x 0 is carried by a set K in 
the family .ft'(A). Since K is a compact set in A, then by using the Krein-Milman 
Theorem we obtain a net {p>o}a of discrete measures converging to /i 0 in A. Using 
hypothesis (HI), we can easily extend each discrete initial measure /Xq to a discrete 
measure p Q in X\ we just apply (HI) to each point in the support of /Xq. By con¬ 
struction, each p a is a tight measure carried by n^ 1 /! n U, which by hypothesis (H3) 
is a compact set. This implies that {pa} a is a uniformly tight net and then Theorem 
2.2 is applied to obtain a subnet converging to some tight measure p, also carried by 
It^K HU, which is, in particular, a Borel set in X. Thus p is a trajectory statistical 
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solution, in the sense of Definition 3.1. The fact that p satisfies the initial condi¬ 
tion, i.e., fl to p = p 0 , follows easily from the uniqueness of the limits in M(X, tight), 
guaranteed by Theorem 2.1. 

The proof of the case when /i 0 is not carried by any set K e .ft'(X), can be reduced 
to the previous case by using the hypothesis that /io is a tight measure and thus, 
in particular, inner regular with respect to the family M'(X) by hypothesis (H2). 
The idea, then, consists in decomposing p 0 as a sum of Borel measures, each being 
carried by a set in &'(X). The previous case can be applied to each of these measures 
(normalized to probability measures), yielding a countable family of W-trajectory 
statistical solutions. Our desired measure is therefore obtained as a (weighted) sum 
of these particular measures. 

In the previous discussion, we skipped some technical details concerning the restric¬ 
tion of the approximating measures to convenient compact subsets. If we assumed 
that our underlying phase space was completely regular, the proof could be made a bit 
simpler, as these restrictions would no longer be necessary since in completely regular 
Hausdorff spaces the weak-star semi-continuity topology coincides with the weak-star 
topology (see Lemma 2.1). But again, looking for a higher degree of generality, we 
assume only that our phase space X is a Hausdorff space. 


Proof of Theorem 3.1. Let us first suppose that /io is carried by a compact set K 
in .h ; (A"). Using the Krein-Milman Theorem ([64, Theorem 3.23]) we obtain a net 
{Hq } a of discrete measures in V(K) such that /Xq w ^* po\ k- Since each /xq is a discrete 
measure, there exist J a G IN, 0 < 9f < 1, and xxq ■ £ K , such that 


3 = 1 




with ~ U f° r every a. 

From the hypothesis (HI) it follows that, for each Uqj, there exists uf E U such 
that n to uf = Uqj. Consider the measure p a defined on X by 


Pa — ^ Oj S uf ■ 

3 = 1 

Note that p a belongs to V(X, tight) and is carried by n DU, which is a compact 
set by the hypothesis (H3). Thus, {pa} a is clearly a uniformly tight net. By Theorem 
2.2, there is a measure p in V(X , tight) such that, by passing to a subnet if necessary, 


p a W -^* p in X. (35) 

Moreover, using Lemma 2.1, we find that p is carried by n^ 1 /! fl hi. 

Consider a bounded and upper semicontinuous function ip : X —> R. Then tp o n fo 
is a bounded function on X. Moreover, since n to is continuous, p o n io is also an 
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upper semicontinuous function on X. Applying a change of variables as in (8) and 

qjj 

using the convergence p a —*■ p together with Lemma 2.1, we obtain that 


limsup / tp(u) dU to p a (u) = limsup / (ip o n to )(w) dp Q (u) 
a J X a J X 


< / (pn f0 )(li) dp(u) = / <p(u) dHtoP(u). 


1 X 


IX 


This means, using Lemma 2.1 once again, that IL 0 p Q Yl to p in X. Further¬ 
more, taking the restrictions of these measures to the compact K , we also have that 
II toPalx U ’^* n to p\ K . On the other hand, we have by construction that 

n to p a \ K = p« w -^* p 0 \ K . 


Adding this to the fact that U to p\K, po \k G V(K, tight), we obtain, by the uniqueness 
of the limit guaranteed by Theorem 2.1, that n to p|^ = /i 0 |a'- Then, since LI to p and 
/i 0 are carried by K we find that Ll fo p = p Q . 

We have thus obtained a measure p G V(X, tight) carried by the compact and 
hence Borel set IT^ 1 A DU C U, with the initial condition Ll fo p = p 0 , which means 
that we have proved the existence of a trajectory statistical solution satisfying the 
initial condition in the case that /i 0 is carried by a set K G 8!(X). 

Now let us consider the case when /i 0 is not carried by any set K G 8!(X). Since 
Po is a tight Borel probability measure on X, we have, by the hypothesis (H2), that 
p 0 is inner regular with respect to the family X(X). Thus, there exists a sequence 
{K n } n of sets in M'(X) such that 

p 0 (K n+ 1 ) > p 0 (K n ) > 0, and p 0 (X\K n ) < - , Vn G IN. (36) 

n 

Moreover, we may assume that K n C K n+1 , for all n G IN. 

Let D\ = K\ and D n = K n \K n _ { , for every n > 2. Note that 

p 0 [j D^j = p 0 (x\ \J K^j < p 0 (X\K n ) < i 


for all n G IN. Thus, taking the limit as n —> oo above, we obtain that po is carried 
by u, Dj. Then, for every A G since the sets Dj, j G IN, are pairwise disjoint, 
we have 


p 0 (A) = p 0 




^p 0 (An Dj). 

3 =1 


ho £ho(^W 

3 


So we may decompose p 0 as 
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where p J 0 is the Borel probability measure dehned as 

Note that each /i q is well-dehned, since /lo(-Di) = p 0 (Ki) > 0 and 


ho (Dj) = Ho(Kj) - > 0 , Vj > 2. 

Also, since each p J (j is carried by the set Kj e .^'(X), using the first part of the proof, 
for each j e IN we obtain a tight Borel probability measure pj carried by Tl^Kj D U 
and such that II to Pj = p J 0 . Let then p be the Borel probability measure dehned by 


Observe that 



P = J2po( D j)Pr 
j 

£ho(^(n to Kj n u) — y ' p,p(Dj) = l, 

3 3 


where the first and second equalities follow from the fact that pj is carried by II ^ K 3 0 
U. Thus, p is carried by (J. B^ 1 K :I P\ U, which is a Borel set in X and is contained in 
U. The fact that Lh 0 p = /i 0 is also easily verified. 

It only remains to show that p is a tight measure. In order to prove so, consider a 
Borel set A E IB x an d e > 0. Let n G IN be such that 1/n < e/2. Since pj is a tight 
measure, there exists, for each 1 < j < n, a compact set /C” C A such that 

pM\k?) < 

Let /C n = Ui an d note that 

OO 

p(A\ic‘) = 

3 =1 

n oo 

< J] MO (Dj) 

j= 1 j=n+l 

< | + MoA\A'„). 

Thus, according to (36) and the choice of n, it follows that p(A\IC n ) < e. Since each 
KA is a compact set in X, this proves that p is tight. □ 


Remark 3.6. Notice that given an initial tight Borel probability measure po on 
X, if /i 0 is carried by a set K e if'(X), then the trajectory statistical solution p 
with n to p = po obtained in the proof of Theorem 3.1 is carried by the Borel set 
Il^K ft U.. On the other hand, if po is not carried by any set K e JT(X), then given 
any sequence of sets K n in .h'(X), n G IN, such that po(X \ K n ) —> 0, as n —> oo, 
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a trajectory statistical solution p with II to p = po can be constructed such that it is 
carried by the Borel set U D {f] n I\.f^K n ). 

The hypotheses in Theorem 3.1 are needed for the solvability of the initial value 
problem for trajectory statistical solutions for arbitrary initial tight Borel probability 
measures on X. It may happen, however, that the conditions are met only for a 
subset of tight Borel probability measures on X, or, similarly, that the continuity of 
the solutions only hold in a coarser topology than that in the phase space. In this 
regard, the same argument of the proof of Theorem 3.1 yields the following result, 
which we state without further details. 

Theorem 3.2. Let X be a Hausdorff space and let I be a real interval closed and 
bounded on the left with left end point to- Let U C X be a subset of X = C\ oc (I,X) 
and let X 0 be a Borel subset of X. Suppose that 
(HI’) U to U D X 0 . 

Suppose, also, that there exists a family &'(X 0 ) of compact subsets of X 0 such that 
(H2’) Every tight Borel probability measure po on X which is carried by X 0 is inner 
regular with respect to the family X(X 0 ) in the sense of (6); 

(H3’) For every K G ^'(X 0 ), the subset n^ 1 K C\U is compact in X. 

Then, for any tight Borel probability measure po on X which is carried by X 0 , there 
exists a U-trajectory statistical solution p on I such that n to p = po- 

Remark 3.7. The extension of the result to the framework described in Theorem 3.2 
is motivated by the original framework of Vishik and Fursikov for the Navier-Stokes 
equations (see e.g. [76]), in which, instead of using the weak continuity in L 2 of the 
weak solutions, one uses the continuity in sufficiently large negative powers of the 
Stokes operator. In this case, one can take, for instance, X = V' = D{A~ 1 ^ 2 ) and 
X 0 = H w , where H and V are as in Section 4.1. 

3.3. Existence of Statistical Solutions in Phase Space. We start with the fol¬ 
lowing lemma. 

Lemma 3.1. Let X be a Hausdorff space, L C R be an interval, and X = C\ oc (L,X). 
Let p be a Borel probability measure on X. Suppose that Z is a Borel subset of X 
and that p is carried by a Borel subset V of the space Z defined in (26). Then, the 
projected measure p t = H t p is carried by Z , for almost every t G /. 

Proof. Consider the characteristic function \z of Z in X. Since V is such that, for 
all u G V, we have u{t) G Z, for almost every t G /, then \z{u{t )) = 1, for almost 
every t G /. Thus, 

[ [ xz(u{t))dtdp(u)= f [ ldtdp(u) = (t 2 ~ti)p(V), (37) 

Jv Jt 1 JV Jt 1 

for every ti,t 2 G /, ti < t 2 . Notice that the map (t,u) h -> Xz(u(t)) is the composition 
of the function xz with the evaluation operator U(t,u ) = u(t) given in (30). Since Z 
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is a Borel set in X , then xz is a Borel function on X. Moreover, since the evaluation 
operator U is continuous and hence Borel from / x X into X, it follows that the 
composition map (t,u) H > xz(u(t)) is also a real-valued Borel function on / x X. 
Then, we apply Tonelli’s Theorem ([25, Theorem III.11.14]) to the left-hand side of 
(37) to obtain that 


-t 2 


Xz(u(t)) d p{u) d t = (t 2 - ti)p(V), 


for all ti,t 2 G /, t\ < t 2 . Using that p is a probability measure carried by V, we find 
that 

[ Xz{u(t)) d p(u) = 1, 

Jv 

for almost every t G I. Using again that p is carried by V, we see that 


Pt{Z) = [ Xz(v)p t (v) = [ Xz(u(t))p(u) = [ xz(u(t)) d p(u) = 1, 

Jx Jx Jv 

for almost every t G I. Therefore, since p t is a probability measure, pt{X \ Z) — 0, 
and hence p t is carried by Z. for almost every 1 G /. □ 


Now we prove that the family of measures obtained as the projections of a trajectory 
statistical solution at each t G / on X is a statistical solution, in the sense of Definition 


3.2. 


Theorem 3.3. Let X and Z be Hausdorff spaces and Y be a topological vector space 
such that Z C X C Y' v *, with continuous injections, where Y' denotes the dual space 
of Y and Yf, is the dual space Y' endowed with the weak-star topology. Consider 
an interval I C R and a subset U C X, where X = C\ oc (I,X). Let p be a U- 
trajectory statistical solution and let V be a Borel subset of X such that V C U and 
p(V) = 1. Assume that 03 z C 03x/ that U C X\, where X\ is defined in (27); and 
that F : / x Z —>• Y' is an (£/ <S> 53y /)-measurable function such that (28) holds 

and ut = F(t,u ) in the weak sense (25). Assume, moreover, that 


1 1 —y 

for every v G Y. Then, 


\(F(t,u(t)),v) Y ',y\ dp(ti) 6 l! oc (/), 


(38) 


<1 >(u(t))dp(u) = / <&{u{t'))dp{u) + / (F(s,u(s)),^'(u(s))) Y ',ydp(u)ds, (39) 


t> Jv 


for all t,t' G / and for all cylindrical test function $ inY'. Moreover, the function 


it-)- / < p(u{t )) d p(u) (40) 

Jv 

is continuous on I for every p G C^(X). In particular, the family of projections 
{pt}tei, where p t = hi t p, is a statistical solution in phase space of the equation u t = 
F(t,u), over the interval I, in the sense of Definition 3.2. 
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Proof. Since V C U C X\, it follows by the definition of X\ that, for every u e V, 
u{t) e Z , for almost every tel. Since C 23 x, we have, in particular, that Z is 
a Borel subset of X. Therefore, it follows from Lemma 3.1 that pt = Htp is carried 
by Z , for almost every tel. This means that p is carried by Uf 1 Z, for almost every 
tel. Hence, the integrals in (38) and (39) with respect to p in V, with integrands 
containing the mapping u e V ha F(t,u(t)), are well-defined almost everywhere in /. 
Now consider a cylindrical test function $ : Y' —> R given by 

$(w) = <j)({u, vi) y >,y, ••■,(«, v k ) Y 'y), Vw e Y', 

where 0 is a continuously-differentiable real-valued function on H k with compact 
support, k e N, and Vi ,..., v k e Y. 

For every u 6W, since U is included in the space 3d defined in (21), the function 
t (—)■ (u(t ), Vj)y' ,y is absolutely continuous on / and, for almost every tel, 

^(u(t),Vj) Y ',Y = (F(t, u(t)),Vj) Y i jY , Vj = l,...,k. 

Thus, 

d ^ ^ 

— <l>(u(t)) = dj(f)((u(t), Ui)y/,y, . . . , {u(t),V k ) Y ’,Y)-^{u(t),V j ) Y ’,Y 

3 =1 
k 

= Y V ')y',Y, • • • , (u{t),V k ) Y ', Y ){F(t, u(t)),Vj) Y > iY 

3 = 1 

= (F(t,u{t)),&(u(t))) Y > tY , (41) 


where is the differential of $ in Y', given in (29). 

Let us show that, for every u eU, the mapping t H > <E>(w(i)) is absolutely continuous 
on /. Since each djf is bounded in R fc , there exists M > 0 such that ||V0(x)|| < M, 
for every x G R fc , where || • || denotes the norm in R fc . Then, given any finite sequence 
of pairwise disjoint subintervals {(tj, Sj)}jLi in /, we obtain, from the Mean Value 
Theorem, 

N 

Y |0((«(Sj)^l)w,F, • • • , (u(Sj), V k ) Y ' }Y ) - <t>({u{tj),v i)y/,y, • • • , («(tj), V k ) Y ',v)\ 

3 =1 

N 

< M E ||((li(Sj),Ul)y',Y, . . . , (li(s J -),Ufc)y/ ) y)-((lt(t J -),U 1 )Y' ) Y, • • • , V k ) Y ', Y ) || • 
1=1 

Thus, the absolute continuity of the mapping 1 H > <f>(w(t)) follows by using that each 
mapping 1 1 —> (u(t),Vj) Y ' yY is absolutely continuous, for j — 1,..., k. 
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Therefore, from (41) we obtain that 

<f>(w(f)) = $(«(£')) + [ { F (s,u(s)),^\u(s)))Y',yds, 
Jt' 


for every t,t' G / and every u EU. 
Consider the function 


( 42 ) 


H : I xX -> R 

(t,u) H* $(w(f)). 


Denote by t the continuous injection of A" into and let U : / x X — > X be the 
function defined in (30). Then H can be written as the composition $ o i o U. And 
since $, t and U are continuous functions, so is H. Similarly, we obtain that the 
mapping (t,u) H > <f>'(w(t)) is continuous on / x X. 

Using the function G : l x X — > Y' defined in (32), with W — Y’ and w = 0, we 
write 


(G(t, u), &(u(t))) Y >,Y 


(F(t, u(t)), $'(u(t)))rr,y, if u{t) G z, 

0, if u[t) G X \ Z. 


Since F is (£/ <g) *Bz, 23y/)-measurable by hypothesis, we have from Proposition 2.1 
that G is (£/ <g) 05^, Q3y')-measurable. Therefore, using also the continuity of the 
function p : Y’ x Y —y R given by 


p(u,v) = (u,v) Y ’,Y, 

it follows that the mapping (t,u) H y ( G(t,u ), $ , (u(t)))y/ i y is (£/ <g) 23 ,y ) - measurable. 
Moreover, recall that Q'(u(t)) is of the form (see (29)) 

k 

$'(«(*)) = V ^)y',Y, •••,(«, Ufe)y/,y)Uj, 

1=1 

where k G IN, V\,.. ., G Y, and 0 is a continuously differentiable real-valued func¬ 
tion on R fc with compact support. Thus, 

k 

{F(t, u(t)), &(u(t))) Y ',Y = dj(j)((u, V\) Y ',Y , •••,(«, V k ) Y ', y){ F ( t, u(t)), Vj) Y ',Y■ 

1=1 

Since the derivatives dj(J) are uniformly bounded, we have that 

k 

\{F(t,u(t)),&(u(t))) Y >,Y\ < IWX*)Wv',v}| > 

l=i 

for suitable constants C\,... ,Ck- Using hypothesis (38) and Tonelli’s Theorem on 
nonnegative measurable functions ([25, Theorem III. 11.14]), we then find that the 
map 


(t,u) 1 y (F(t,u(t)),&(u(t))) Y ',Y 
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is integrable on / x X, with respect to the product of the Lebesgue measure and the 
measure p. 

Thus, from (42), and considering that V C U is a Borel subset of X, we obtain 

r r r rt 


4>(w(t)) dp(u) = / <f>(w(f')) dp(w) + / / (G(s,u),&(u(s))) y >,y dsdp(u), (43) 


Jv Jv JvJt' 

for every f, t' G I. Applying Fubini’s Theorem ([25, Theorem 111.11.9]) to the second 
term on the right hand side of (43) we find that 


<&(u(t))dp(u) = / $(u(t'))dp(u) + / / (G(s, u), &(u(s)))Y',Ydp(u)ds 


t< Jv 


= / §{u(i!))dp{u) + / / (F(s > u(s)) ) ^ , (u(s)))yf iy dp(u)ds, (44) 


Jv Jt' Jv 

for all f, t' G /. This proves the mean equation (39). 

Now consider a function (p G C^(X) and let us prove that the function defined in 
(40) is continuous on /. Given i G /, consider a sequence {t n } n in / such that t n —$■ t. 
Since every u G V is continuous from / into X , it follows that 


p(u(t n )) ip(u(t)), Mu G V. 


Since p is in particular a bounded function on X, we obtain, from the Lebesgue 
Dominated Convergence Theorem, that 

/ <p(u{t n )) d p(u) -t / <p(tt(t)) dp(it). 

Jv Jv 

Thus, the function defined in (40) is continuous on I. 

Concerning the family {pt}tei of the projections p t = I l t p, tel, it is immediate 
to see that condition (i) of Definition 3.2 follows from the continuity of (40), that 
condition (ii) of Definition 3.2 follows from hypothesis (38), and that condition (iii) of 
Definition 3.2 follows from the mean equation (39). Therefore, {pt}tci is a statistical 
solution in phase space. □ 


The next result provides a solution for the Initial Value Problem 3.2. Given an 
initial measure /i 0 on X , we use Theorem 3.1 to obtain a trajectory statistical solution, 
which is then projected at each time t to yield a family of measures on X. Thanks 
to Theorem 3.3, this family of projections is a statistical solution. 


Theorem 3.4. Let X and Z be Hausdorff spaces and Y be a topological vector space 
such that Z C X C Yh* , with continuous injections, where Y' denotes the dual space 
of Y and Yf* is the dual space Y' endowed with the weak-star topology. Let I C R be 
an interval closed and bounded on the left with left end point to and let U be a subset 
of X = C\ oc (I,X). Suppose that U satisfies the hypotheses (HI), (H2), and (H3) of 
Theorem 3.1, and assume, moreover, that the following conditions hold 
(H4) C 
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(H5) U C X\, where X\ is defined in (27), and F : I xZ —>Y' is an (£/(E>©Z) 93y')- 
measurable function such that (28) holds, and u t = F(t, u ) in the weak sense 
(25); 

(H 6 ) there exists a function 7 : / x X x Y —» R such that, for every v G Y, the 
function ( t,u ) 1 —> 7( t,u,v ) is (£/ <E) Q 5 x) -measurable and 



I (F(s,u(s)),v) Y ',y 


ds < y(t,u(t 0 ),v), 


Vt el, Vue U. 


Then, for any tight Borel probability measure po in X such that 


(45) 


l(t,u 0 ,v) dfi a (u„) < oo, 


lx 


for a.e. t e /, and for all v eY, 


(46) 


there exists a projected statistical solution {p t } te j of u t = F(t,u), over the interval I, 
associated with aU-trajectory statistical solution, such that p to = p 0 . 


Proof. Since U is a subset of X satisfying all the hypotheses of Theorem 3.1 and p 0 
is a tight Borel probability measure on X, it follows from that theorem that there 
exists a 77-trajectory statistical solution p on / such that n to p = po- 

Since p is a 77-trajectory statistical solution, there exists a subset V of U such 
that V G VSx and p(V) = 1. As in the proof of Theorem 3.3, it follows from the 
hypotheses (H4) and (H5) that the integrand on the left hand side of (45) is integrable 
with respect to the product of the Lebesgue measure on / and the measure p on X. 
Then, we take the integral in (45) with respect to p and apply Tonelli’s Theorem ([25, 
Theorem 111.11.14]) to find that 



|(F(s,u(s)),n)y/ i y| dp(w)ds < / y(t,u(t 0 ),v) dp(u) 


t 0 Jv 


= / y(t,u 0 ,v) dU to p(u 0 ) = / y(t,u 0 ,v) dp 0 (u 0 ). 


’X 


lx 


Thus from (46) we obtain that 


[ \( F (t,u(t)),v) Y ',Y\ d p(u) e L\ oc (I). (47) 

Jv 

This proves condition (38), so that, together with (H5), F is an (£/ ® 03^, 93y')- 
measurable function satisfying all the hypotheses of Theorem 3.3. All the other 
hypotheses on U and on p of Theorem 3.3 are already assumed to be satisfied. Then, 
it follows from that theorem that the family {p t }tei, with p t = HtP, is a statistical 
solution satisfying p to = hb 0 p = p 0 . □ 


Similarly, we may consider an extension of Theorem 3.4 with the set of hypotheses 
(HI’), (H2 7 ), and (H3 7 ) described in Theorem 3.2. In this case, we have the following 
result, which we state without further comments. 


ABSTRACT FRAMEWORK FOR THE THEORY OF STATISTICAL SOLUTIONS 


35 


Theorem 3.5. Consider the framework of Theorem 3.4 and letX 0 be a Borel subset of 
X. Suppose thatU satisfies, instead, the conditions (HI’), (H2’), and (H3’) described 
in Theorem 3.2. Assume, moreover, that the conditions (H4), (H5) and (H6) of 
Theorem 3-4 hold. Then, for any tight Borel probability measure po in X which is 
carried by X 0 and satisfies (46), there exists a projected statistical solution {pt}tei, 
associated with a U-trajectory statistical solution, such that p to = /io- 

Remark 3.8. Concerning the hypothesis (H6) and the condition (46) of Theorem 
3.4, we will, in fact, have, in many applications, that Y is a Banach space, and we 
will obtain the stronger condition that there exists an (£/® 93x)-measurable function 
7 : / x X —> R such that 

[ ||F(s, w(s))||y' ds < 7 ft,uft 0 )), Vf £ I, \/u G IA. (48) 

Jto 

In this case, the existence of a projected statistical solution for the initial value 
problem will hold for any initial Borel probability measure /i 0 such that 


/ y(t, Mo) dp 0 (uo) < oo, for a.e. t G I, (49) 

Jx 

In the following result, we obtain a mean energy inequality for trajectory statistical 
solutions, provided the individual solutions satisfy a corresponding energy inequality. 


Proposition 3.1. Let X, Z, Y as in the statement of Theorem 3.3. Consider U C X\, 
where X\ is defined in (27). Let p be a U-trajectory statistical solution and let V be 
a Borel subset of X = C\ oc (I,X) such that V CU and pfV) = 1. Suppose that there 
exist functions a : L x X —$■ R and (3 : I x Z ^ R which are (£/ ® 93 x, 93 r.) and 
(£/ <g) 53z, 93 r) measurable, respectively, and satisfy the following conditions 

(i) (t, u ) i— y aft, uft)) belongs to P(Jx V, Ax p), for every compact subset J C I; 

(ii) ft, u ) i— y (3ft, uft)) belongs to L l { J x V, A x p ), for every compact subset J C I; 

(iii) For p-almost every u G V it follows that 

^a(t,uft)) +P(t,u(t)) < 0, (50) 

in the sense of distributions on I, i.e., 


ip'(s)a(s,u(s))ds + / ip(s)/3(s,u(s))ds < 0, 


for all nonnegative test functions p G 


(51) 


I Jv 


p'(s)a(s, u(s))dp(u)ds + / / ip(s)j3(s, u(s))dp(u)ds < 0, (52) 


i Jv 


Then, 
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for all nonnegative test functions p G Cf°(I). In terms of the projected statistical 
solution in phase space {pt}t&i, where p t = II t p, we can write (52) as 


I JV 


p'(s)a(s, u)dp s (u)ds + / / p(s)/3(s,u)dp s (u)ds < 0. 


(53) 


/ JV 


Proof. The proof follows by integrating (51) with respect to p on V and then applying 
Fubini’s Theorem ([25, Theorem III.11.9]) by using hypotheses (i) and (ii). □ 


The motivation for considering the passage from the inequality (50), valid for in¬ 
dividual weak solutions, to the mean inequality (52) comes from the Navier-Stokes 
equations (see Section 4.1) and other similar equations from fluid flows. It also ap¬ 
pears in different types of equations, such as the nonlinear wave equation considered 
in Section 4.3. In some situations, however, such as in the case of the Reaction- 
Diffusion equation considered in Section 4.2, the individual weak solutions satisfies in 
fact an energy-type equality, and of course this equality is similarly passed on to the 
statistical solutions. In the case of equality, the test functions are allowed to assume 
negative values. We state this result as follows, omitting the proof since it follows 
along the same lines as that of Proposition 3.1. 


Proposition 3.2. Under the hypotheses of Proposition 3.1, if the equality holds in 
(50), then the equality holds in (52) and (53), for any test function p G as 

well. 


4. Applications 

In this section we apply the general framework developed previously to obtain the 
existence of statistical solutions for particular examples of evolution equations. 

4.1. Incompressible Navier-Stokes Equations. The Navier-Stokes equations are 
a commonly used model in the study of turbulent flows. In their three-dimensional 
form for incompressible Newtonian fluids, these equations are written as 

du 

— — uAu + (u • V)u + Vp = f, (54) 

V • u = 0, (55) 

where u = (iq, m 2 , W 3 ) is the velocity held, p is the kinematic pressure, f represents a 
given body force applied to the fluid, and v is the parameter of kinematic viscosity. 
We consider u, p and f as functions of a space variable x and a time variable t, with 
x varying in a set 0 C R 3 and t varying in an interval / C R. 

For a physical formulation of the equations we refer the reader to the books by 
Landau and Lifshitz [50] and Batchelor [ 6 ]. For the mathematical approach, see 
Ladyzhenskaya [48], Temam [69, 70] and Constantin and Foias [20]. 

Our main concern in this section is to apply the abstract theory developed in the 
previous sections to prove the existence of statistical solutions to the Navier-Stokes 
equations (54)-(55). For more specific discussions on the notion of statistical solutions 
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for the Navier-Stokes equations, we refer the reader to [27, 75, 76, 33] and also to the 
more recent paper [36]. Except for the fact that we allow for a slightly less regular 
external forcing term, the existence result that we obtain in this section has been 
already proved in [36], a work that has been in fact our inspiration. Nevertheless, 
we want to show that our abstract framework does apply to this case, with a much 
simpler proof than that presented in [36]. 

For the sake of simplicity, we assume periodic boundary conditions. In this case 
we consider a periodic domain given by Q = nf =1 (0, Li), where L* > 0, for i = 
1, 2, 3. This means we are assuming the flow is periodic with period L { in each spatial 
direction x*. We also consider the averages of the flow and of the forcing term to be 


/ u(x, t) dx = 0 , / f(x, f) dx 

Jn Jn 


= 0 . 


Let R 3 ) denote the space of infinitely differentiable and fl-periodic functions 

u. We define the set of periodic test functions with vanishing average and divergence 
free as 


V°° := ^u£C r (f!,R : 


V • u = 0 and / u(x) dx = 0 

Jn 


(56) 


Let H be the closure of V°° in L 2 (h2,R 3 ) and let V be the closure of V°° in 
H 1 ( h2,R 3 ). The inner product and norm in H are defined, respectively, by 


where u ■ v 


(u, v) = / u ■ v dx and |u| = \J (u, u), 
■hi 

hli=\ u i v i■ hi the space V, these are defined as 


((u, v)) = (Vu, Vv) = / Vu : Vv dx and ||u|| = \J ((u, u)), 

Jn 

where it is understood that Vu = (dui/dxj)f ■ =1 and that Vu : Vv is the componen¬ 
twise product between Vu and Vv. We also consider the space H endowed with its 
weak topology and denote it by H w . 

Clearly, V is a subset of H . Thus, by identifying H with its dual space H\ we 
obtain the following continuous inclusions 

V C H = H' cV'. 


Since the injection H ‘-A V' is a continuous linear mapping, we also have 

H w > C 

with continuous injection, where Vh denotes the space V' endowed with the weak-star 
topology. 

Let A be the Stokes operator, defined as A = —PA, where P : L 2 (h2,M 3 ) — > H is 
the Leray-Helmholtz projection, i.e., the orthogonal projector in L 2 (fi,M 3 ) onto the 
subspace of divergence-free vector fields. We denote by D(A) the domain of A, which 
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is defined as the set of functions u G V such that Alu G 77. In the periodic case with 
zero average, we have 

Au = -Au, Vu G D(A) = V n 77 2 (fi, M 3 ) 

and A is a positive self-adjoint linear operator with compact inverse, so that it has 
a nondecreasing sequence {A*}^ of positive eigenvalues counted according to their 
multiplicity, associated with an orthonormal basis {wj} ie ]N in 77 . Furthermore, the 
Poincare inequality holds, i.e., for all u G V, 

Ai|u| 2 < H 2 , (57) 

where Ai > 0 is the first eigenvalue of the Stokes operator. 

The Stokes operator A extends to an operator defined on V, with values in V', and 
so that 

||Au||y/ = sup {A\l,v) V iy = SUp ((u, v)) = ||u||, 

||v||=l ||v||=l 

which implies in particular that A-.V—tV'is continuous. 

We denote by P k : 77 —* V the Galerkin projector onto the space spanned by the 
eigenfunctions associated with the first k eigenvalues of the Stokes operator, i.e. 

k 

Pk U = y>, Wj)wj, Vu G H. 

i= 1 

Since {wj} ie ]N is an orthonormal basis in 77, we obtain 

k k 

Hull 2 = (0PfcU,P fc u)) = l( u > w i)l 2 (( w i> w i)) = l(u> W J’)I% < ^fc|u| 2 , 

3 = 1 j =1 

for every u G 77, from which it follows that Pk is continuous from 77 to V. 

Let B : V x V —> V' be the bilinear operator 

7?(u, v) = P[(u ■ V)v], Vu, v G V, 

which satisfies the inequality (see [20, 33]) 

||P(u, v)H < c|u| 1 / 4 ||u|| 3 / 4 |v| 1 / 4 ||v|| 3 / 4 , Vu, v G V, (58) 

where c is a universal constant. By using this inequality, it is not difficult to see that 
B : V x V —* V' is also a continuous operator. 

The natural space for the solutions of the Navier-Stokes equations is the space 
Ci oc (7, 77 w ) of continuous functions from an interval 7 C R to 77 w , endowed with 
the compact-open topology. This function space is the space of weakly continuous 
functions from 7 to 77. 

The Navier-Stokes equations can be written in the functional form 

u t + uAu + 77(u, u) = f in V', 


(59) 
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ill the sense of Bochner (see Remark 2.6), which in this case is equivalent to the weak 
formulation 

^(u, v) + z/((u, v)) + ((u • V)u, v) = (f, v), 

in the distribution sense on /, for any v 6 h 

The notion of solution that is considered here is the well-known Leray-Hopf weak 
solution, which is defined below. 

Definition 4.1. Let I be an interval in R and f G Lf oc (I,V'). We say that u is a 

Leray-Hopf weak solution of the Navier-Stokes equations (54) -(55) on I if 

(i) U e H) n p oc (/, v) n C loc (/, h,) ; 

(ii) s t u 6 V'); 

(iii) u satisfies the functional formulation of the Navier-Stokes equations 

Uf + vAu. + £>(u, u) = f; (60) 

(iv) u satisfies the energy inequality in the sense that for almost all t' G I and for 
all t E I with t > t', 

^|u(t)| 2 + z 'J ||u(s)|| 2 ds < \\ n W)\ 2 + J (f(s),u(s)) V / )V ds; (61) 

(v) If I is closed and bounded on the left, with left end point to, then the solution 
is strongly continuous in H at t 0 from the right, i.e., u(t) —>■ u(t 0 ) in H as 
t^-tf. 

The set of allowed times t' in (61) can be characterized as the points of strong 
continuity from the right of u in H. In particular, condition (v) implies that t' = to 
is allowed in that case. 

The Leray-Hopf weak solutions of the Navier-Stokes equations also satisfy a strength¬ 
ened form of the energy inequality (61) of Definition 4.1. The proof of this strength¬ 
ened energy inequality has been given in [33] for external forces f in Ll c (I,H). This 
inequality is also valid if f belongs to the larger space L 2 oc (J, V '), as we state below: 

Proposition 4.1. Let T > 0 and f G L 2 (0,T;H / ). Consider a nonnegative, non¬ 
decreasing and continuously-differentiable real-valued function if : [0, oo) —$■ R with 
bounded derivative. If u is a Leray-Hopf weak solution of the Navier-Stokes equations 
on [0, T], then 

^WH*)! 2 )) < 2 V' / (|u(t)| : - : )[(f(t),u {t))vy ~ H|u(i)|| 2 ] 
in the sense of distributions on [0, T]. 

The idea of the proof is to first obtain such inequality for mollifications of if and 
|u(t)| 2 , and then to pass to the limit with respect to the parameters of each mollifi¬ 
cation, in a suitable order. For the detailed proof, see [53]. 
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Given R > 0, we denote by Bh(R) the closed ball centered at the origin and with 
radius R in H. The corresponding closed ball endowed with the weak topology is 
denoted by B H (R) W . We define the following sets of Leray-Hopf weak solutions: 


Ui = {u E Ci oc (I, H w ) : u is a Leray-Hopf weak solution on /}, (62) 

Uj(R) = {u E Cio C (7, Bh(R)w) : u is a Leray-Hopf weak solution on /}, (63) 

U\ — { u E Ci oc (J, H w ) : u is a Leray-Hopf weak solution on /}, (64) 

U\(R) = {u G Ci oc (J, B h (R) w ) : u is a Leray-Hopf weak solution on /}, (65) 

where I denotes the interior of the interval I. 

From (61), one obtains the classical estimates 

|u(t)| 2 < KOI 2 + ^||f||i a (66) 
[ ||u(s)|| 2 ds <-|u(0| 2 +4ll f lli 2 (t',t;E') ; ( 67 ) 

Jt' v v 

valid for any uGh/, for any time t' G / allowed in (61) and for any t G / with t > t'. 
Moreover, using (58), one also has that 

( pt \ 3/4 5/4 

jf \Ms)\\$6s) < ^i|u( 0| 2 + ( 68 ) 


where c is a universal constant and D(t',t) is a nondimensional function which de¬ 
pends on the variables t',t and also on the parameters u, Ai, A2, and f through the 
nondimensional quantities u\i\t — t'\ and (A/ 1 /K 2 )||f|| L 2 {t’,t-y)- 

The a priori estimates ( 66 )-( 68 ) allow us to prove that U\(R) is a compact and 
metrizable space, in the same way as it was done in [36, Proposition 2.2], Furthermore, 
one can show that li\[K) is the closure of the space Ui(R) with respect to the topology 
of C(J,if w ). 

The existence of a Leray-Hopf weak solution on a given interval I C R is obtained 
using the estimates ( 66 ), (67), and ( 68 ). This proof is a classical result and can be 
found in many well-known texts [20, 48, 52, 70]. We state it below for completeness. 

Theorem 4.1. Let I C 1R be an interval closed and bounded on the left with left end 
point t 0 and let f G L 2 oc (J, V'). Then, given u 0 G H, there exists at least one weak 
solution u eUi of (54)-(55) in the sense of Definition f.l satisfying n to u = u 0 . 

From now on, we assume that / C R is an interval closed and bounded on the left, 
with left end point to- Under this assumption, the energy inequality (61) is valid for 
t' = t 0 . 

Consider a compact subinterval J C I. Then given u G 77/ such that u(to) E Bh(R) 
for some R > 0, it follows, from ( 66 ) with t' = to, that there exists R > R such that 
u(t) G B h {R), for every t G J. Thus, the restriction of u to J belongs to Uj(R). 
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In order to prove the existence of a trajectory statistical solution for the Navier- 
Stokes equations satisfying a given initial data, we shall apply Theorem 3.1 by con¬ 
sidering X as the space H w and the general set U as the set of weak solutions Ui. We 
now show that Ui satisfies the hypotheses of Theorem 3.1. 

First, note that hypothesis (HI) is a direct consequence of Theorem 4.1. Also, 
defining as the family of (strongly) compact sets in H , hypothesis (H2) follows 

from the fact that H is a separable Banach space (see Sections 2.2 and 2.3). The 
only hypothesis that needs more work is the remaining hypothesis (H3), which is the 
subject of the next proposition. 

Proposition 4.2. Let I C R, be an interval closed and bounded on the left with left 
end point t 0 and let K be a strongly compact set in H. Then HUj is a compact 

set in X = Ci oc (J, H w ). 

Proof. Let u G n Uj and let R > 0 be sufficiently large so that K C B H (R). 

Consider an increasing sequence {J n } n of compact subintervals of / with left end 
point at t 0 and such that / = (J n J n . Since n (o u G K C B H (R ) and u 6 Zf;, it 
follows, from the estimate ( 66 ) with t' = t 0 , that there exists a sequence {R n }n of 
positive real numbers such that n^u G Uj n (R n ), for every n. Thus, 

npKnUjd^n-j'UjSRu). (69) 

n 

Since each Uj n (R n ) is a metrizable space, (69) implies that n HUi is also metriz- 
able. Therefore, it suffices to show that P\Ui is sequentially compact. 

Let {ufc}fc be a sequence in n HUi. As in the classical proof of existence of 
weak solutions (Theorem 4.1), using the a priori estimates ( 66 )-( 68 ) on each compact 
interval J n and applying a diagonalization method, we obtain a subsequence {u^/}*./ 
and a function u such that 

Life' y u in Ci oc (J, H w ), (70) 

as k' —> oo. Moreover, this limit function u is a weak solution on the interior of /, 
i.e. u G Li\ (the condition of strong continuity at to, item (v) of Definition 4.1, is not 
guaranteed at this point). From (70) we obtain in particular that 

Ufc'(to) ->■ u (t 0 ) in H w . (71) 

On the other hand, since K is a compact set in H, there exists a further subsequence, 
which we still denote by {u*,/}*/, and an element u 0 G K such that 

u fc /(t 0 ) -t u 0 in H. (72) 

From (71) and (72) it follows that u(to) = Uo, which implies that u G More¬ 

over, we obtain that {u fc /(t 0 )}fc' also converges to u(t 0 ) in the strong topology of 
H. 

This allow us to prove that u verifies in addition the last condition of Definition 4.1. 
Indeed, since each u fc / belongs to Uj, they satisfy in particular the energy inequality 
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(61) at t' — t 0 - Considering the liminf as k' — > oo in this inequality we obtain, by 
using the strong convergence of {u fc /(f 0 )} fc / to u(f 0 ) in H , that 

\\u(t )\ 2 + v f ||u(s)|| 2 ds < i|u(t 0 )| ;J + f (f(s),u(s))v/ >v ds. 

Z JtQ Z Jt 0 

Then, by taking the lim sup as t — * tf above, we obtain 

lirnsup |u(£)| 2 < |u(t 0 )| 2 . (73) 

Since u e Ci oc (J, H w ), we also have that 

|u(t 0 )| 2 < liminf |u(£)| 2 . (74) 

t^t+ 

Now (73) and (74) imply that u(t) converges in norm to u(to) as t —> tf. Since u(t) 
also converges weakly to u(f 0 ) as t —> , we deduce that 

lim u(£) = u(£q) in H. 

This means that u G hl io 1 A" D Ui, which completes the proof of compactness of this 
set. □ 

Now we are able to prove the existence of a solution for the corresponding Initial 
Value Problem 3.1 associated with the Navier-Stokes equations. 

Theorem 4.2. Let I C R be an interval closed and bounded on the left with left 
end point t 0 and let Uj be the set of Leray-Hopf weak solutions of the Navier-Stokes 
equations on I. Then, given a Borel probability measure p 0 on H, there exists a Ul¬ 
tra] ectory statistical solution p on C\ oc (I , H w ) satisfying the initial condition II to p = 
ho- 

Proof. Our intention is to apply Theorem 3.1 to the set Uj , with X = H w and 8 !(H W ) 
as the family of (strongly) compact sets in H. First of all, since H is a separable 
Banach space, the Borel sets in H and H w coincide (see Section 2.3). This implies 
that any Borel probability measure p,Q on H is also a Borel probability measure on 
H w , and vice-versa, so that we can refer indistinguishably to probability measures 
on either H or H w . Moreover, since H is a Polish space, it follows that any Borel 
probability measure on H is tight in the sense of being inner regular with respect 
to the family of compact subsets of H ([1, Theorem 12.7]). Thus, 8 !(H W ) satisfies 
the hypothesis (H2) of Theorem 3.1. From Theorem 4.1, it follows that the set Uj 
satisfies hypothesis (HI) of Theorem 3.1. From Proposition 4.2, we also obtain that 
Uj satisfies hypothesis (H3) of Theorem 3.1. Thus, Uj verifies all the hypotheses 
of Theorem 3.1. Therefore, we apply Theorem 3.1 to deduce that there exists a 
77/-trajectory statistical solution p with n to p = p Q . □ 
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Finally, using Theorem 4.2 and the strengthened energy inequality from Proposition 
4.1, we obtain a solution for the corresponding Initial Value Problem 3.2 associated 
with the Navier-Stokes equations. More specifically, we prove the existence of a 
projected statistical solution of the Navier-Stokes equations, in the sense of Definition 
3.3, associated with a Ui -trajectory statistical solution and satisfying a given initial 
data. 

Theorem 4.3. Let I C R be an interval closed and bounded on the left with left 
end point t 0 and let Uj be the set of Leray-Hopf weak solutions of the Navier-Stokes 
equations on I. Consider a Borel probability measure p 0 on H satisfying 


L 


u | 2 d/ro(u) < oo. 


(75) 


Then there exists a projected statistical solution {pt}tei of the Navier-Stokes equations 


(54)-(55) ? associated with a Ui-trajectory statistical solution, such that 

(i) The initial condition p to = po holds; 

(ii) The function 



is continuous on I, for every bounded and weakly-continuous real-valued func¬ 
tion ip on H, and is measurable on I, for every bounded and continuous real¬ 
valued function ip on H. 

(iii) For any cylindrical test function <h in V', it follows that 




(f(s) — vA\\ — B(u, u), $'(u))v/',y dp s (u)ds, (77) 


+ / / f s - 


for all t, t' e I. 


(iv) The mean strengthened energy inequality 



(78) 


is satisfied in the distribution sense on I, for every nonnegative, nondecreasing 
and continuously-differentiable real-valued function 0 with bounded derivative. 

(v) At the initial time, the limit 



(79) 


holds for every function 0 as in (iv) . 
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Proof. We have seen in the proof of Theorem 4.2 that the set of Leray-Hopf weak 


solutions Hi satisfies all the hypotheses (HI), (H2), and (H3) of Theorem 3.1. Now 
let F : / x V —* V' be the function defined by 


F (t, u) = f (t) — vA\i — B( u, u). 


As previously mentioned, the linear operator A : V — > V' and the bilinear operator 
B : V x V —> V' are continuous. This implies that the mapping u ha — i/Au — R(u, u) 
is also continuous from V into V'. In particular, the mapping (t, u) ha — vAu— R(u, u) 
is (£/ <E) 23y, 23 y,)-measurable. Further, since f G L^ oc (I,V'), we obtain that F is a 
(£/ <8) 23 y, 23 v 7 )- measurable function. From the functional equation (59) and the fact 
that any weak solution u 6 Wj satisfies tq G L 4 / 3 (/;F'), it follows that (28) holds. 
From the condition (iii) of the Definition 4.1 of a Leray-Hopf weak solution, the 
validity of u t = F(t, u(£)) in the sense of Bochner implies that (25) holds in the weak 
sense (see Remark 2.6). Hence, hypothesis (H5) of Theorem 3.4 holds. 

The a priori estimate ( 68 ) means that there exists a function 7 : / x H w -A R such 
that 



/ ||F(s, u(s))||v' ds < 7(t,u(t 0 )), Vt G /, Vu G Uj, 
Jt 0 

which is clearly (£/ ® 23^ w )-measurable and, thanks to (75), with 



Thus, hypothesis (H 6 ) and condition (46) hold (see Remark 3.8). 

Also, as mentioned before, we know that V C H w C IA*, with all the injections 
being continuous. Since H is a separable Banach space, then 23 h w = 23 h- Moreover, 
since V is a Polish space, then 23 y C 23# = 18# w (see Section 2.3), showing that 
hypothesis (H4) is also verified. 

Then, applying Theorem 3.4 with X = H W} Z = Y = V, U = Ui , F and 7 as above, 
we obtain the existence of a projected statistical solution {p t } teI associated with a 
Uj -trajectory statistical solution p and such that p to = /io- This means that {p t }tei 
satisfies (i), (iii), and the first part of (ii), concerning bounded and weakly-continuous 
functions ip on H. 

Let us prove the second part of property (ii), concerning strongly continuous func¬ 
tions. Consider a bounded and continuous real-valued function ip on H. Let P m , 
m G IN, be the Galerkin projectors. Then, for every m G IN, the function p o P m is 
bounded and continuous on H w . Let V C U be a Borel subset such that p(V) = 1. 
From the first part of (ii), it follows that the function 



is continuous on /, for every m G IN. Then, since the function (76) is the pointwise 
(in t) limit of these functions asm A 00 , it follows that (76) is measurable on I. This 
proves the second part of (ii). 
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For the proof of (iv), consider the functions a : / x H w —> R and j3 : / x V —> R 
dehned respectively by 

oi(t,u(t)) = ip (|u(f)| 2 ), (80) 

and 

P(t,u(t)) = -2^(\u(t)\ 2 )[(f(t),u(t))vy - H|u(*)|| 2 ], (81) 

for every u G 77/ and t G I. Using the Galerkin projector as above, we see that a and 
/? are the pointwise limit of continuous functions, hence they are measurable maps as 
required in Proposition 3.1. Using the estimates (66) and (67) with t' = to, which is 
allowed for functions in Uj, and using (75), we obtain that a G L°°(J x V, A x p) and 
/3 G L 1 ( J x V, A x p), for every compact subset J C I, where A denotes the Lebesgue 
measure on /. From Proposition 4.1, the functions a and /3 satisfy 

-^o;(f, u(t)) + P(t, u(t)) < 0, Vu gWj, 

in the sense of distributions in /. Property (iv) then follows by applying Proposition 
3.1 with X = H w , Y = V,U=Ui and with the functions a and f3 dehned in (80)-(81). 

ft only remains to prove property (v). Note that for every function i/j as in (iv), 
we may write 

^(|u(t)| 2 )<^(0 )+^(O|u(t)| 2 , 

for some 0 < £ < |u(t) . Then, using the boundedness of -0' and the a priori estimate 
(66) with t' — to, we find that 

^(|u(t)| 2 )<Go + C' 1 |u(t 0 )| 2 , 

for suitable constants C 0 ,Ci > 0. Hence, from (75), it follows that , 0(|u(t); J ) is 
bounded by a p-integrable function in V which does not depend on t. Furthermore, 
since every u G Uj is strongly continuous at t 0 and 0 is continuous, we have that 

0(|u(f)| :J ) -»• 0(|u(f o )| 2 ), 

p-almost everywhere, as t —> t^. Therefore, (79) follows from the Lebesgue Dominated 
Convergence Theorem. □ 

4.2. Reaction-Diffusion Equation. In this section, we consider the following reaction- 
diffusion-type equation 

Ou 

— (x, t) — aA«(x, t) — f(t,u(x, t)) + p(x, t), x G fl C R", tel, (82) 

subject to the boundary condition 

tt(x,t)| xe an = 0, Vi G I, (83) 

where u is the unknown variable, a is a positive constant, / is the reaction function 
and g is the external force. Moreover, / C R is an arbitrary interval and fl C R" is 
a bounded and open subset which is assumed to be smooth. 
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We follow the same framework and notations from [18, Section XV.3], but in order 
to simplify the presentation we consider only a scalar equation instead of a system of 
equations. 

Consider the spaces H = L 2 {fl) and V = H/(fl) with respective norms | • \h and 
|| • \\ v , given by 

\ v \h = ( |u(x )| 2 dx, Wv G H, 

Jn 

and 

||v|[v = i |Vu(x)| 2 dx, WveV. 

Jn 

Also, consider V' = H^^fl), the dual of H/(fl), with duality product (•, -)v',v- Then, 
identifying H with its dual space H ', we have 

V C H = H' C V\ 

with continuous inclusions and, in particular, H w Vf* with continuous injection. 
We also consider the space ff), for r > 0, and its dual H~ r (fT), with H/(fl) C H C 
H~ r (fl), for every r > 0. For r > n/2 — n/p and p > 2, we have H r 0 (Q) C L P (Q). We 
denote the duality product between L P (Q) and L q (Q), 1 < p, q < oo, 1/p + 1/q = 1, 
simply by (•,•), which includes the inner product of H. 

We assume that g G L 2 oc (I , V ') and that / is a function in C(R x R, R) satisfying 
the following estimates, for every v G R and s G R: 


V\v\ p - < f(s,v)v, 

(84) 

|/(s,u)|^T < C 2 {\v\ P + 1), 

(85) 


where 77 > 0, p > 2 and C\, C 2 G R are constants. In [18], the function g is also 
assumed to be translation bounded in the space T 2 oc (J, V') 1 but we do not make this 
assumption since we do not need uniform estimates for arbitrarily large times. 

As in [18], it follows by using condition (85) that if r > max{l,n(l /2 — 1/p)} and 
u G Lf oc (J, LP(ft)) n L 2 oc (/, V), then d t u G L£ c (/, H~ r (Q)), for 1/p + 1/q = 1. This 
implies that the evolution equation (82) can be considered in the distribution sense 
on /, with values in H~ r (Q). 

We then have the following definition of a weak solution for problem (82)-(83). 

Definition 4.2. A weak solution of (82)-(83) is a function u = u(x,t) on fix I such 
that u G L p oc (I, L p (fl)) D L 2 oc (J, V) and u satisfies (82) in the distribution sense on 
I, with values in H~ r (fl). 

Given R > 0 , let Bh(R) be the closed ball centered at the origin and of radius R 
in H. Consider the following sets of weak solutions: 

Uj = {u G C loc (I, H) | u is a weak solution of (82)-(83) on /}, ( 86 ) 

Ui(R) = {uG C\ oc (J, B h (R)) I u is a weak solution of (82)-(83) on /}. (87) 
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The proof of existence of individual weak solutions for the corresponding initial 
value problem of (82)-(83) can be found in [18]. We state it below for completeness. 

Theorem 4.4. Consider an interval I C R bounded and closed on the left with left 
end point t 0 . Let g G Lf oc (I, V') and let f G C(R x R, R) be a function satisfying 
conditions (84) and (85). Then, given uq G H, there exists a weak solution u of 
problem (82)-(83) such thatu G Lf oc (J, L p (Ll))C\L 2 oc (I, V)nL“ c (J, H) andu(t 0 ) = u 0 . 

The following proposition presents some additional properties satisfied by every 
weak solution of (82)-(83) in the sense of Definition 4.2. The proof is given in [18, 
Proposition XV.3.1]. 

Proposition 4.3. Let u G Lf oc (J, L p (fl)) nL 2 oc (I, V) be a weak solution of (82)-(83). 
Then 

(i) u G C\ oc (L, H); 

(ii) the function \u(s)\ 2 H is absolutely continuous on every compact subinterval 
J C I and satisfies the following energy equality 

^K*)l h + a|K*)llv + (f(t,u(t)),u(t)) = (< g{t),u{t)) v >y , (88) 

for almost every t G /. 

Note that, in ( 88 ), we abuse notation by denoting as f(t,u(t )) the mapping x G 
fl hg f(t,u(x,t)), which, as a consequence of (85), belongs to L q {kl) for almost every 
t G / and for all u G Ui. 

Now we prove that the set of weak solutions Ui satisfies the hypotheses of The¬ 
orem 3.1. We first observe that Theorem 4.4 implies that Ht 0 ldj = H , so that 1Ai 
satisfies hypothesis (HI) of Theorem 3.1, with X = H. For the remaining hypotheses 
of Theorem 3.1, we actually prove the stronger property that (H3) holds for every 
compact subset of X. This is given in the following proposition. Hence, in this case, 
we can take 8!(X) to be the family of all compact subsets of X, so that, in particular, 
hypothesis (H2) is trivially true. 

Proposition 4.4. Let I C R be an interval closed and bounded on the left with left 
end point t 0 and let K be a compact subset of H. Then nD Ui is compact in 
Ci oc (I,H). 

Proof. Since X = C\ oc (I, H) is a metrizable space, it suffices to show that CUi 

is sequentially compact. Consider then a sequence { Uj} :j in n ^. 1 K n Uj. Since K 
is compact there exists uq G K such that, by taking a subsequence if necessary, 
iij(t o) —>■ uq in H. This implies in particular that the sequence {tq(lo)} is bounded 

in H. 
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Using condition (84) on the energy equality ( 88 ) for each Uj and integrating from 
t 0 to t, we obtain that 

\uj{t )\h ~ \uj(t 0 )\ 2 h + a [ \\ Uj (s)r v ds + 2rj [ \uj(s)\ p LV ds 

Jto Jto 

<-/ ||s'(s)||y, ds + 2Ci\t - t 0 ||ft|, (89) 
a Jto 

where | • \u> denotes the norm in L P (Q), |0| is the Lebesgue measure of hi, and Cj 
is as in (84). Consider a sequence {J n }n of compact subintervals of I such that 
/ = (J n J n . Then, from the estimate (89) and the boundedness of the sequence 
{uj(t 0 )} in H, it follows that, for each n, {uj}j is a bounded sequence in L 2 (J n , V) D 
L p (J n , L P (Q)) D L°°(J n , H). Using the same arguments as in [18, Theorem XV.3.1] 
and a diagonalization process, we obtain a weak solution u of problem (82)-(83) on 
/ such that, modulo a subsequence, Uj —> u in C\ oc (I,H). In particular, it follows 
that Uj(t ) —> u(t ) in H , for every t G I. Thus, w(to) = Uq E K and we conclude that 
u G n^K nh/, as required. □ 

The existence of a trajectory statistical solution with respect to a given initial data 
now follows by a simple application of Theorem 3.1 for X = H and Ui as the set of 
weak solutions of (82)-(83) over a given interval / C R closed and bounded on the 
left. Recall that since H is a Polish space then every Borel probability measure on 
H is tight. We then have the following result. 

Theorem 4.5. Let I C R be an interval closed and bounded on the left with left end 
point t 0 and letUi be the set of weak solutions of problem (82)-(83) on I. If p 0 is a 
Borel probability measure on H then there exists a Uj-trajectory statistical solution p 
on Cioc(C H) such that II to p = po- 

Next we obtain the existence of statistical solutions with respect to a given initial 
measure. 


Theorem 4.6. Let I C R be an interval closed and bounded on the left with left end 
point to and let p 0 be a Borel probability measure on H satisfying 

/ \u\ 2 h dpo(u) < oo. (90) 

J H 

Then there exists a projected statistical solution {p t }tei of (82)-(83), associated with 
a Ui-trajectory statistical solution, such that 

(i) The initial condition p to = po holds; 

(ii) The function 

t f ip{u) dp t {u) (91) 

J H 

is continuous on I for every bounded and continuous real-valued function <p 
on H; 
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(iii) For any cylindrical test function $ in H r (f2) ? with r > max{l,n(l/2 — l/p)}, 
it follows that 


$(m) d p t (u) = / <E>(w) d p t '(u) 


i H 


IH 


+ 



f Jh 


(aAu - f(s , u) + g(s), &(u)) H - r ( n),jj5(n) dp s (-u)ds, (92) 


for all t, f G I; 

(iv) For every nonnegative, nondecreasing and continuously-differentiable real-valued 
function 0 with bounded derivative, the function 


1 1 ^ 


ip(\u\ 2 H ) d p t {u) 


’ H 


is absolutely continuous on I, and the following mean strengthened energy 
equality holds in the distribution sense on I, 


d 

dt 


1 H 


0 (M h) d Pt{u) = 2 fj\\u\ 2 H )[{g(t),u)v,v ~ a\\u\\ 2 v - ( f(t,u),u )] dp t {u). 


IH 


Proof. The proof follows by arguments similar to the ones used in Theorem 4.3. We 
apply Theorem 3.4 and Proposition 3.2 with X = Ft, U = Uj, Z = Hq( Q) fl L P (Q), 
Y = Hq(Q), and the functions F : / x Z —)■ Y 1 , at : I x H —)■ 1R and (3 : I x V —» R 
defined as 

F(t, u ) = aAuit) - fit, u(t )) + g{t), 


a{t,u) =fj(\u\ 2 H ), 

and 

P(t,u) = -2'ip\\u\ 2 H )[{g(t),u) v , y - a\\u\\y - (. f(t,u),u )]. 

□ 

4.3. Nonlinear wave equation. In this section, we apply the abstract framework 
to prove the existence of statistical solutions of a nonlinear hyperbolic-type equation 
which appears within the theory of Relativistic Quantum Mechanics. We follow the 
framework presented in [52, Chap. 1, Sec.l]. 

Let O C R n be a bounded open set with smooth boundary, denoted by dfl, and let 
/ C R be an arbitrary interval. Consider the equation 

(Fii 

- Au+\u\ r u = f, (93) 

where the real-valued function u = u(x.,t) is the unknown variable, r is a positive 
constant and / = /(x, t ) is a given function, with x G 12 and t G /. 

We endow equation (93) with the following boundary condition: 

“(x,*)lxesn = 0, Vie/. 


(94) 
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In order to obtain a functional setting for problem (93)-(94), we introduce the space 

v = ^(n)nf/(fi), 

where p = r + 2. The space V turns into a Banach space when endowed with the 
v ’ 


norm || • 11 T >, dehned by 


\v 


— IMI/rt + \ v \ lp , Vf G V, 


where || • ||#i and | • \u> denote the usual norms in the spaces Hq(D,) and L P (Q), 
respectively. 

The dual space of V is the space V' = /7 _1 (h2) + L 9 (il), where l/p+ 1/q = 1. The 
duality product between V and V' is denoted by (•, -)y, y. 

Also, we consider the space L 2 (h2) endowed with its usual norm and inner product, 
which are denoted respectively by I - I r2 and (-, -)r 2 - Moreover, we assume that f is a 
function in L 2 oc (/,L 2 (ft)). 

We rewrite equation (93) in the following equivalent form: 

du 


at ” = 0 ' 

^-A«+i«r u = i 


(95) 


We denote the nonlinear term of the second equation in (95) by the function b : 
V V' given by 

b{u) = \u\ r u, Vu G V. 

Then, considering U = (u, v ) and the linear operator A dehned by 


AU = 

the system (95) becomes 


0 -I 
-A 0 


U = 


0 -I 
-A 0 


—v 
—A u 


dU 

d t 


+ AU + N(U) = G, 


(96) 


where N(U) and G are the vectors 

N(U) = 


0 

b(u] 


G = 


From (94) we obtain the following boundary condition for (96): 

£/(x,t) | xe an = 0, Vt e I. (97) 

Now we define V — V x L 2 (h2), with the norm 

ll^llv = yJ\\A 2 v+\v\v, V£/=MgV. 

When endowed with its corresponding weak topology, the space V is denoted by 
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We characterize the dual of V as the space V' = L 2 x V' (see Remark 4.1), with 
the duality product between h = (/, g ) G V' and U — (u, v) G V as 

(h, U)v,v = (/, v)l* + (g, u)v',v- 

With this representation, the usual norm for an element h = (/, g ) in the dual space 
V' can also be written as 

llfcllv = sj\f\b + llsll \r 

We now give the dehnition of a weak solution of problem (96)-(97). 

Definition 4.3. Let I C R be an interval and let f G L 2 oc (I, L 2 (h2)). We say that 
U = U(t) = (u(t),v(t)) is a weak solution of problem (96)-(97) on I if the following 
conditions are satisfied: 

(i) UeL^(I,V); 

(ii) CeCi O0 (/,Cw); 

(iii) U satisfies 

^ +AU + NiU) = G in V', (98) 

at 

in the sense of distributions on I; 

(iv) For almost every t’ G I, U satisfies the following energy inequality 

E{U{t)) < E(U(t')) + [ ( G{s ), U(s)) V ' tV ds, (99) 

■it’ 

for every t G / with t > t!, where G = (/, 0) and 

E(U) = E(u,v) = l -\\uf Hl + ^\u\ p LP + ^\v\ 2 L 2 ] (100) 

(v) If I is closed and bounded on the left, with left end point to, then U is strongly 
continuous at t 0 from the right, i.e. U(t ) —$■ U(to) inV as t —$■ t£. 

The set of times f for which (99) is valid can be characterized as the points of strong 
continuity from the right of U, and they form a set of total measure in I. 

Remark 4.1. A more natural and usual representation for the dual of V = V x L 2 (Ll) 
is simply V' x L 2 (fl), preserving the order of the variables in the product space. 
However, in order to be able to write the equation for the evolution of the information 
f v $(U) d pt(U) in the general form (33) (see (108)), we switched the representation 
to V' = L 2 ( Q) x V. For the same reason, we introduced the term G = (/, 0). 

For any R > 0, let B V (R) denote the closed ball of radius R in V. We define the 
following sets of weak solutions of problem (96)-(97): 

Ui = {U G C\o C (I, Vw) : U is a weak solution of (96)-(97) on /}, (101) 

Uj(R) = {U G C(I, B V (R) w ) : U is a weak solution of (96)-(97) on /}. (102) 



52 


ANNE C. BRONZI, CECILIA F. MONDAINI, AND RICARDO M. S. ROSA 


Next we state an existence theorem of individual weak solutions for the initial value 
problem associated with the system (96)-(97). The proof is given in [52, Theorem 1.1, 
Chap. 1, Sec. 1]. Although the regularity conditions (ii), (iv) and (v) of Definition 
4.3 are not explicitly written in the statement of the theorem in this reference, they 
are obtained along the lines of their proof. 

Theorem 4.7. Let I C K be an interval closed and bounded on the left with left end 
point t 0 and let f G L 2 oc (J, L 2 (D)). Then, given U 0 G V, there exists at least one weak 
solution U EUj of (96)-(97), in the sense of Definition f.3, satisfying n to U = U 0 . 

Consider now / C R an interval closed and bounded on the left with left end point 
t 0 . In this case, item (v) of Definition 4.3 implies that the energy inequality (99) is 
valid for t' — t 0 . 

Let U = (u,v) G Uj such that U(t 0 ) G B V (R ) , for some R > 0. From the energy 
inequality (99) with t' = to, it follows that 

E(U(t)) < R + ^r f |/(s)|| 2 ds+l f \ v ( s )\l 2 ds, (103) 

z Jt 0 z Jt 0 

for every t G I, which also yields 

hv(t )\ 2 L 2 < R+ 1 [ \f(s )\ 2 L2 ds + | f \v(s )\ 2 L 2 ds. (104) 

Z Z J to Z j to 

Then, given a compact subinterval J C /, by applying Gronwall’s inequality in (104) 
we obtain that \v(-)\ 2 l2 is uniformly bounded on J. From the estimate (103), it then 
follows that there exists R> R such that U(t) G B V (R), for every t G J. Thus, the 
restriction of U to J belongs to Uj(R). 

We shall now prove that the set of weak solutions Ui satisfies the hypotheses (HI), 
(H2), and (H3) of Theorem 3.1. 

Theorem 4.7 shows, in an equivalent form, that n tQ Ui = V m . Thus, the set Ui 
satisfies hypothesis (HI) with X = V w . 

Now define 


8!(V W ) = {K C H w | K is a (strongly) compact set in V}. 

Since V is a separable Banach space, we obtain that every Borel probability measure 
/i 0 on V w is tight with respect to the family &'(V W ) (see Sections 2.2 and 2.3). Then, 
considering X = V w , it follows that the family •ft'(Vw) satisfies hypotheses (H2) of 
Theorem 3.1. The next proposition shows that lAj also satisfies hypothesis (H3). 

Proposition 4.5. Let I C R, be an interval closed and bounded on the left with left 
end point t 0 and let K be a set in Then n^ 1 A" D IA] is a compact set in 

X = Ci oc (I, V w ). 
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Proof. Let R > 0 be such that K C By(R) and let {J n }n be a sequence of compact 
subsets of / such that 

I = U J n • 

n 

Then, from the energy inequality (99) with t' = to, one obtains that, for every n G IN, 
there exists a positive real number R n > R such that 

n Jn UeU Jn (R n ), vu en^KrUj, 

which implies that 

n^xnu, cf|n A '« A W. 

n 

Since Uj n (R n ) is a subset of C\ oc (I, By(R n ) w ), which is a metrizable space, it follows 
that II^K nW/ is also metrizable. Thus, it is enough to prove that ft Uj is a 

sequentially compact space. 

Consider then a sequence {Uk}k in D Ui. Since Uk(t 0 ) G K and K is a 

compact set in V, there exists Uq G V and a subsequence {kj}j such that 

Ukj(t o) —t Uq in V. (105) 

Following classical arguments used for the existence of weak solutions (see [52, 
Chap. 1, Sec.l]), we obtain a priori estimates that allow us to pass to the limit 
on each compact set J n . Then, using a diagonalization process, we obtain a further 
subsequence (which we still denote by {Uk j }j) and a function U defined on the interval 
/ such that {Uk }j converges to U in Ci oc (I,V w ) and U is a weak solution on the 
interior of I. Thanks to (105) we have at the initial time that [/(to) — Uq e K, so 
that U G fl^iC Then, as in the case of the Navier-Stokes equations (see the proof of 
Proposition 4.2), using the energy inequality and the fact that the convergence (105) 
at the initial time is in the strong topology, we obtain that U is strongly continuous 
at the initial time t 0 , so that U G Ui. Therefore, U G Ll^ 1 /! D Ui, proving that 
P\Uj is compact. □ 

Thus, applying Theorem 3.1 with X = V w and U as the set of weak solutions Ui , 
we obtain the following result on the existence of a trajectory statistical solution for 
the nonlinear wave equation with respect to a given initial data. 

Theorem 4.8. Let I C R be an interval closed and bounded on the left with left end 
point t 0 and letUj be the set of weak solutions of problem (96)-(97) on I. If /i 0 is a 
Borel probability measure on V then there exists a Ui-trajectory statistical solution p 
on Ci oc (I, V w ) such that fI fo p = p 0 . 

In the next result we prove the existence of a statistical solution of problem (96)- 
(97) in the sense of Definition 3.2 with respect to a given initial data. 
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Theorem 4.9. Let I C R be an interval closed and bounded on the left with left end 
point t 0 and let Uj be the set of weak solutions of problem (96)-(97) on I. Consider 
a Borel probability measure /i 0 on V satisfying 


E(U) dfio(U) < oo, 


(106) 


Jv 

with E as defined in (100). Then there exists a projected statistical solution {pt}tei 
of (93), associated with a Hi-trajectory statistical solution, such that 

(i) The initial condition p to = po holds; 

(ii) The function 


t ha [ <p{U) d p t (U) (107) 

Jv 

is continuous on I, for every bounded and weakly-continuous real-valued func¬ 
tion ip on V, and is measurable on I, for every bounded and continuous real¬ 
valued function on V; 

(iii) For any cylindrical test function $ in V', it follows that 


[ <*(£/) d p,(U) = 

Jv 

[ $(U)dp t ,(U)+ f [ (G-AU-N(U),&(U)) V ',vdp 8 (U)ds, (108) 
Jv Jt’ Jv 

for all t, t' G I; 

(iv) The mean strengthened energy inequality 

4 J ip(E(U)) d p,(U) < J V(E(U))(G(t), u(t)) v .,v d Pl (U) (109) 

is satisfied in the distribution sense on I, for every nonnegative, nondecreasing 
and continuously-differentiable real-valued function if with bounded derivative, 
where G = (/, 0); 

(v) At the initial time, the limit 

lim f if{E{U)) d pt(U) = [ f>{E{U)) d/i 0 (U) (110) 

t-ttf Jv Jv 

holds for every function if as in (iv) . 

Proof. The proof follows by arguments similar to the ones used in Theorem 4.3, 
considering X = V w , Z — Y = V and the functions F : I x V —> V', a : I x V w —> R 
and /? : / x V —* R defined respectively as 

F(t, U ) = G(t) - AU(t) - N(U(t)), 

a (t, U) = if(E(U(t))), 


P(t,U) = -if'(E(Um(G(t),U(t )) v , y , 


and 
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for every (t, U ) in the corresponding domains. □ 
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